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INTRODUCTION 

The technical uteclal presented in this report describes the work 
under NASA Grant NSG-1603, Langley Research Center, Haapton, VA for the 

period of Kerch 16, 1981 -> October 30, 1962, Dr. Eugene D. Denaep served 

li 

ea principal Investigator for the grant. 

The work under this grant was to investigate numerical algorithms for 
large-^space structures with particular emphasis on decoupling itethod for 
analysis and design. The luvsstigator and project participants have con- 
sidered numerous aspects of the analysis of large systems ranging from the 
algebraic theory to lambda matrices to lde!^..tiflcatlon algorithms. Previous 
reports have described soma of the problems considered including the theory 
of the decoupling procedure. The material presented in this report is a 
general treatment of the algebraic theory of lambda matrices and application 
of the theory to second-order lambda matrices. Since the finite-element 
analysis of the dynamics of a structure can be characterized by a second- 
order matrix differential equation, the second-order lambda matrix is of 
Importance in analysis and design. 

Although published literature on lambda matrices is adequate in some 
aspects, there are numerous gaps in the general theory and a serious lack 
of computer software for analysis of large structures. As an example, there 
is no software package known to this investigator that computes the latent 
roots and latent veuitors with efficiency comparable to the software available 
for computing eigenvalues and eigenvectors of general matrices. The usual 
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computational approach for handling lambda matrlcaa la to form a atata- 
varlabla matrix from the lambda cuitrix and then uae the exlatlng aoftware 
In computet llbtarles. Since the atate-variable matrix la aparae, the exla- 
tlng software may not be efficient for such problems. 

Chapter 2 of thli report glvep the general theory of lambda matrices. 

This work was carried out to acquire a better understanding of lambda matrices 

* 

as well as fill In some of the gaps In the theory. The relationships between 
eig.iinvalues » latent roots, eigenvectors, latent vectors, elgenprojectors ahd 
latent projectors are presented In this chapter. The decoupling algorithm is 
based on Rlccatl matrices and projectors so this chapter plays an Important 
role in the work that follows. The theory of lambda matrices presented In the 
chapter should also help in developing better aoftware, as wall as reference 
material for other researchers. 

A lambda matrix of the general form A(l)"AQX”+A^X*'’*^-f , . .+A^ Is said 
to be regular If det(AQ)^0. Chapter 2 Is restricted to the regular lambda 
matrix as the dynamic of a structure can generally be defined In the regular 
form. Lambda matrices which are not regular have been considered in Chapter 
3. This type of lambda matrix will sometimes result from the analysis of 
Che dynamics of a system. A fifth order lambda matrix Which is not regular 
was encountered in some related work on the space shuttle. The theory In 
Chapter 3 complements the development in Chapter 2 and the two chapters 
cover the general theory of lambda matrice regardless of regularity. 

The material in Chapter 4 of the report is not directly related Co the 
large-space structure problem, but is given for completeness as well as for 
the researcher working with low-order as well as low dimensioned system. The 
procedure for computing latent projectors or matrix residues given in that 
chapter would not be preferred to the eigenvector method given In Chapter 2 
and 3 unless the worker has the inverse lambda matrix, tA(X)l 
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The flnlte-eleqienit: model fchet fchie inveatlgetor hee conaldMred Is cher- 
ecterlzed by e second-order lambda matrix. Chapter 5 presents soma theorama 
on second-order lambda matrices and details on the elgenptoj actors and latent 
projectors of such lambda matrices. The material presented In the chapter 6 

can be used as a basis for ccsaputlng the damping matrix although an In-depth 
Investigation reveals that the computational burden for this method of coan* 
putlng the daR^lng matrix will be high. Work will continue on determining the 
damping matrix by use of the mathematics of this chapter. 

The final chapter of the report describes a decoupled optimal control 
algorithm tha appears to be reasonably efficient for large space structures. 

The method discussed a decoupling procedure which can be used with optimal 
control theory to damp selected modes In a struccure., A general software 
package for determining feedback control for a structure could be developed 
from the theory. Investigations made to date do not Indicate a limitation on 
the procedure other than the normal onef^ encountered In large systems. The 
mode spill-over problem is eliminated by the method. 

The material presented In this report required the effort of several 
Individuals. The contributions made by Dr. Graham Goodwin, Jesus Leyva-Samos 
and G.Z. Jeon have been essential to the Investigations and the material pre- 
sented in this report. Dr. Graham Goodwin was on leave from the University 
of Newcastle, New South Wales, Australia and was a Visiting Professor for the 
six month period January- June 1982 at the University of Houston. The contri- 
bution to Chapter 3 is hereby acknowledged as well as the time given to research 
personnel for discussion of various problems. 

The Principal Investigator is grateful to NASA Langley Research Center 
for the support during the grant period as well as the encouragemsnt and 
suggestions of Dr. Garnett Homer, technical monitor of the grant. 


CHAPTER 2 


THE ALGEBRAIC THEORY OF UTEMT VECTORS AND PROJECTORS IN LAMBDA MATRICES 

Multivariable syatema and controls are often fomulated in terms of n-th 
order matrix differential equations which Rive rise to lambda matrices of the 
form A(X) ■ A^X”+Aj^ . »+A^. This chapter describes the algebraic 

theory of latent roots | latent vectors » and latent projectors end gives the 
relationships to eigenvalues, eigenvectors and elgenprojectors of the com- 
panion form matrlXt The theory presented in the chapter Is general In that 
distinct as well as repeated eigenvalues or latent roots are treated. The 
work Is limited to regular lambda matrices, det (A^) 0« 

The chapter consists of three sections. Section 1 Is Introductory In 
scope and presents the mathematical forms of lambda matrices, and companion 
matrices as well as the connection between elgenprojectors and latent proctors. 
Section 2 gives the algebraic theory of lament roots, eigenvalues, latent 
vectors and eigenvectors. It Is shown that the latent roots obtained from 
the lambda matrices are equivalent to the eigenvalues are computed from the 
companion matrix. This result la well known and presents no new material. 

The relationships that exists between the left and right eigenvectors and 
the left and right latent vectors Is then developed and It Is shown that the 
latent vectors are subvectors of the eigenvectors. The algebraic theory of 
the elgenprojectors and the latent projectors la developed In Section 3. 

It Is shown chat the latent projectors are submatrlces of the elgenprojectors. 
Section 3 Is a general treatise on elgenprojectors and latent projectors where 
distinct as well as repeated eigenvalues have been considered. 
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1. Introduction 

Multiyarinbl* aystoau auch aa a flnlta-altaant an dal of vibrating a true- 

turaa, control ayatana» and larga ac«la ayataaa can ba formulatad In tantf 

(( 

of aacond or higher ordar natrlx dl'^farantlal aquatlona^ Although auch 
syataaa can ba raformulatad In otata-varlabla for:f,| It oay ba aura officiant, 
froai a nuanrlcal vlavpolnt, to analyze tha ayatan ualng the hlghar-ordar 
dlffarantlal aquatlona. 

To IJ^lustratOi aaouaa that tha dyneaULca of a ayatan can ba eharaetariaod 
by tha nonhonoganaouo antrlx dlffarantlal aquation 

A, . ... + A « • f Ct) (1.1) 

° dt“ ^ dt*^ ^ " 

whom x<t)eR* and f(t)eR*, fron which it followa that A(a) given 

by 

A(a) - AqO^+Aj^o®"^ + ... ... + A^ (1.2) 


will raault whan the Laplace tranafora la taken of (1.1). If a ■ X than 
(1.2) becomea a lapbda aatrlx Cor aatrlx polynomial). If A^ ■ I, than 
(1.1) can be defined In atate-varlable form and the coi^anlon matrix A^ 
can ba vrlttan aa 


^ " 


0 

0 



I 0 

0 I 

a - a dee 

a a baa 

a a a a a 

"Vi “Vi • • ' 



(1.3) 


J) 
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which hM •IgcnvaluM end «lg«nv«cto7i will ba a anXan Aatrlx 

with no aiganvaluaa and no aiganvacto]^a. 

Whan tha partial fraction axpanaion of [A(X)]~^ is takan tha following 
ganaral fon ia obtainad 


[A(x>r^ 


i 


i»l 


*r‘»i p, . 

A-0 (X-X^)^^ 


(1.4) 


whare r la tha nuobar of diatinct latant roota given by dot A(X) ■ 0, la 
tha aultipllcity of a root «id la aqual to tha nullity of A(X^) for tha 
latant root X^. Tha natricaa P^^ cc“ »1U ba callad latant projactora 
(or aatrlx raalduaa, [1]). 

Similarly, tha c(;^mpanlon fora A (X) ia givan by [IX~A 1 which haa aigan- 

6 6 

valuaa X^^^ givan by dat[A^(X)] m 0. Tha partial fraction axpanalon of A~^(X) 
can ba wrlttan aa 


[A^(X)1 


-1 


r 

- I 

1-1 


i-0 


‘i.a 


<X-X^) 


SA-l 


with P^ denoted aa elgenprojactora or aatrlx reaiduea [2], 

The purpose of this chapter is to fonaulata tha algebraic theory of lambda 
aatrlces and the relationship of latant roots, latent vectors and latant 
projectors to tha elgenvaluaa, algenvectora and aigenprojactora of the coa- 
panion fora. The chain rule for latent projectors and elganprojectors for 
tha repeated latant root or eigenvalues will be given. 

This work follows the lines of the earlier work of Lancaster [3], 
Lancaster and Webber [4], and that of Dennis, Traub and Weber [5] on lambda 
nutrlces and aatrlx polynoalala. The reader should refer to Zadeh and Dasoer, 
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[IJ imd CuXlan [6j for tucarlal on utrjbc roiiduao and projactora aa wall 
aa to tha papar by Danaan and Laystla-Raaoa {2] for aoaa aatarlal on aigan- 
projactora. It ia aaauaad that tha raadar la familiar with tha fvuridaiaaRtal! 
of llnaar algabrat An axcallant aourca of ganaral aatarlal on linaar algabra 
ia Gantaaehar, (7]i Lambda iutrleaa ara diacuaaad by Lancaatar in [3], 
and [8]. 

2. Lanbda Hatrlcaa and Coapanlon Foraa 
A ayataa of dlffarantial a<{uationa rapraaantlng a phyaical nodal nay 
ba givan by 

*0 \ * ■ • • V<‘> - 

with initial condltlona x(0) ,x(0) ». . . ,x(0)eR® vhara A^cR®^“, x(t)eR® and 
f(t)eR'*. 'nia axlatanca and unlquanaaa of tha aolutlon to (2.1) la aaauaad. 

It la furthar aaauaad that (2.1) ia Laplaca tranaforaabla with tha Laplaca 
tranafora taking tha ganaral fora 

[AQa“+Aj^a”'-^+...+A^3x(a) - B(i) (2.2) 

whara B(a) containa tha initial condition Infomatlon aa wall aa tha Laplaca 
tranafora of f(t). 

If a 1 in (2.2), than tha rcaultlng aquation can ba conaldarad an a 
lambda matrix aquation. Tha laft hand brackatad tara in (2.2) ia of inter- 
aat in the foaainlng parta of thla papar. Tha aquation in brackata, A(X) 
or A(a), ia fraquantly callad a aatrlx polynoai,al although Lancaatar [3] 
and Dannla, Traub, and Wabbar [5] daflnad A(X) aa a lambda matrix. Thia 
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papav vlll adh«r« to tho tomlnology of L&neutiir , ot. , and xofor to «(X) 
M a laabda aatrlx. 

2.1 Daflnltlon . Glvan a oat of conatant nutritaa and a acalar X 

auch tha'£ XeC» than a laabda aatrlx will ba daflnad aa A(X) ■ AgX'*4>X^X'‘*^4> 

. . Thla natrlac will ba eallad ragular whan dot V 0. For Aq non* 

aingular, than A^ can ba factorad out of A(X) auch that 

A(X) - Ap(lX"+Aj^X““^4-.*,»5'A^) - AqA(X) (2.3) 

2.2 Daflnltlon . Lot A(X) ba aa dafinad in (2.3) and 1st X^oG ha a latant 
root of A(X) whara a latant root X^ la a root of dat A(X) ■ 0. 

If A(X) la ragular . that, thara will ba an latant roota of A(X) which 
ara aqual ta tha alganvaluao of A^, 

2.3 Thaoratft . Lat A(X) ba aa dafinad In (2.1) and lot tha datandnant of 
JT(X) ba glvan by tha ganaral fora 

dat A(X) - CQX"“-f<?j^X"““ V. . .+C^ - 0 (2.4) 


Tha lanbda aatrlx will hava at ooat aa*nulllty Aq latant roota . 
Proof. Lat A(X) ba rapraaantad by 


A(X) - 


.jl(X) 

.jl(X) 

a 

a 


022(X) 

a 

a 


Ota 

a a a 

• • • 
a a a 

a a a 




a..(X) 


(2.5) 
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wh«r« Mch i* * acftlttr polynoalal of «t Mat ordar n, 

•ij "■ ‘ •+*ljn 

vhara a^j^^ ara tha coaffleianta of tha acalar polynoalal with a^^^ ■ 
Tha datarainanfc of (2.5) haa baan dafload by Franklin [9] aa 


dat A(X) ■ ^ a (J . I • • • J_)a« • (X)aM. (X)..«a ^ (X). 

a A j X X ^ ■ X3, 


(2,7) 


Tha auBMtlon axtanda ovar ■! paraueatlona l»2...a with 


a(J|,J2...JJ « aign n 

* X<F<9i» ’ ^ 

Tha coafflciant Cq la obtalnad fron tha landing eoafflclant of dot A(X) aa 
glvan In (2.7) which will ba 

Cq ■ I “idjO *■• *aij^0 ■ ^0 

If dat ■ 0 than A(X) will hava at noat an-1 latant roota. Tha coclflciant 
will ba glvan aa a llnaar cooblnatlon of tha dataralnanta of all tha aub- 
■atrlcaa (b-I)x(b-I) of A^, C2 will ba glvan aa a llnaar conblnatlon of tha 
datamlnanta of all aubaatrlcaa (■-2 )x(b- 2) of A^ and ao on until Aq la tha 
aaro aatrlx. 

2.3,1 Corollary . Thare will ba at laaat nullity A^ latant roota at tha 
origin, i.a. X ■> 0. In tha particular caaa whan dat A^ ■ 0, thara la at 
laaat ona latant root at X - 0. 


ORiaiNW. MCE IS 
’ OF POOR QUALITY 

thm ftooi of fehlo corollary followa dlracrly Uom (2.7) whan tha co- 
afflclanfc In tha iiaMaeion la eonaldarad. 

2.4 DaflniClon . LaC A(X) ba daflnad as In (2,3) with watrlx coafflclonca 
A^CK**** and latent roota of amltlpllclty 'fha nuabar of prlnary right 
or laft latent vactora will ba - nullity A(Aj^). Tha right latent yactora 
will be denoted by y^^ and tha laft by aj^^ where yj^^ and aatiafy 
tha relatione 



J • X»2, ... iq^ 

(2.10) 


i “ X»2,. , . »qj|^ 

(2.11) 

with «d 




Tha prlaary right latent vector y^^ le a aubvector of the linear* in- 
dependent right eigenvector y^^ of the coapanion aatrlx of (1.3) with 






«) 


2„<J) 


Kf 


i'i 


Xn-l-(i) 
Ai y^ 


( 2 . 12 ) 


where y 


(j) 

ci 


aatlaflaa the usual algebraic equation <A^-X^l)y 


<j) 

ei 



* 


Such algenvectora will be referred to ae prlaary aigenvectora in the reat 
of thia paper. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


U 


It follows that ths prlasry right slgsnvsctor Is a function of y^^ 
and X^. Ths ptiaary loft slganvaetor satlsflsa a slallar fora a^; 
with 




. 0 ) . 

*cl 


(XJ“^I+aJxJ"^+ 0.. +-^n-l^*l^^ 
(Xj"^l+A^X"“^+ ... +aJ. 2^*1^^ 


.' i ) 


(2.13) 


whtra la a prlaary latent vector., 

The two foraa of y^^^ and given In (2.12) and (2.13) hold when 
^ j and the oaxlaua nuaber of prlaary right or left latent vectors will 
be a. If > q^ then generalized latent vectors aust be constructed 

. to coapletely define the laabda aatrlx froa Its latent roots and vectors. 
The a^-q^ generalized latent vectors satisfy a chain rule as given In 
Lancaster, [4]. 


2.5 Theorea . Let a(X) be defined as In (2.3), a set of right latent 

f2) (h • ) ]ip<i 

vectors y^ ', y^ ' y^ “* eC fora a right Jordan chain associated 

with the latent root X^ and the jth prlaary latent vector. The chain rule 
la given by 


An)y<^^ + yC2-l) . JL i^>y(^“2> 

'^1 ^ dX ^1 ^21 


dX -"I 

,a-i) 


^ .... 


dX‘ 


( A - l)l 


. 0 

^1 axl 


dX 


C2-1) 


1 * l , 2 ,..., h . 

/ J 


(2.14) 




12 

ORIGINAL RAGE IS 
OF POOR QUALITY 


wh«r« is th« Jt:h prlaary right latent vector end hj le the length of the 
Jordan chain. The vectora for 1 < k 1 hj are the generaliaed right latent 
vectora of the Jth priaury latent vector. 

Proof . The proof of thie theorem is obtained £ton coneideratlon of the chain 
rule for generalized eigenvectore. The chain rule ie 






„( 2 ) 

^ei 


(2.15) 




being the j th primary eigenvector and y^^' for 1 < £ hj its generalized 

eigenvectors associated with the eigenvalue of A^. When (2.15) la expan- 

(k) 

dad the chain rule is obtained where y^ is formed from the first m rows 

* Ck) 

of y^i • 

The generalized right latent vectors y^^^ amd y^^^ will be obtained from 


(k) 


xa )y«> + . 0 

AQAi^yi 

dA(X,) , d^A(X ) 

A(A^)y^ ^ dX ^1 ■ 21 ^2 y^ 


mxl 


or (2.14) in general. 
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Th« chain rule can alao ba utillcad to aodlfy (2*12) for tha relation 
batvaan tha aiganvaetora and tha ganaraliiad latent yactora with 


<k) 

^ci 


*(k) 


X y<k)^<k-l) 

Vi ^i 




1' 

i-0 J * ^ 

lk-i>i] 


(2.16) 


2.6 Thaoram . Let A(^) ba defined aa in (2.3), a aat of left latent yactora 


.(1) -C2) 


.(h.) 


•£ ' » * » » • • • » i 


fom a left Jordan chain aaiociatad with the 
latent root A ^ and tha Jth prlnary latent vector. Tha chain rtila ia given by 

T/^ ^ j2j^Tj 

7^ 


)?a I «- 2 )^ 

A t dX *i ^21 ,,2 *i ^ 


( 2 - 1)1 


d<"-^>A^(X,) 


.( 1 ) 


dX 


(2-1) 


‘'nXl 


i " 1*2, • • • ,h 


‘J 


(2.17) 


where '' ia the jth primary left latent vector and ia tha length of the 


(k^ 

Jordan chain. Tha vactora z^ ' for 1 < k <. h^| ara tha ganarallsad left 
latant vactora of the Jth primary latent vector. 


'^i|/ /y 
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Proof . Th« proof of this thoorom follows dlrsctly from tha gaharsllssd 

Isft slgsnvectors of tha companiton fora, tha ganorallzad laft latant 

(k) (k) 

vactor will ba f oroad from tha laat a rows of z . 

fki 

Tha ganarallzad laft alganvactor z^^' can also ba daflnad fron tha 
latant vactors z^^^ of tha latant roots X^ and tha laobda aatrlx A(X). 
Utilizing (2.13) and tha chain rula for lb followa that tha laft. 

ganarallzad alganyactora satisfy tha ralatlon 



a 

a 

a 

(X Ji+aJx^+Aj) z|‘‘V(2X^I+aJ) 


(X^I+aJ)z^*‘^+z 


(k-1) 

1 



(2.18) 


(k) 

The latant vector z^ ' Is defined only for 1 < k < hj . 

The number of primary latent vectors assoclatad with the latant root 
X^ Is and each of these latent vectors could have a chain of generalized 
latent vactors. Tha structure of a Jordan block *1 could ba given by 
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Xl I 


<2.19) 


1 . 

■uch thAt ^ H. “ jn. » and every Jordan chain haa a very well defined ntmber 

j,l J 1 

of s«nerallzed latent: vectore. The value of h^ la auch that the chain rulee 
for (hj+1) as given In (2.14) and (2.17) are not satisfied. 

An exoaple will now be considered to Illustrate the coaputstlonel pro* 
cedurss. Let A(X) be defined os 

1 ol , f-4.0 2.0 1 fa.o ^1.0 I 

A(X) - M *•' (2.20) 

0 1 -0.2 -3.0J Lo.2 2.6J 


which has latent roots X^ - 1 with amltlpllclty * 1 and X£ 
tlpllclty 82 ■ 3. The right latent vectors are 

V- (y'» y<“ yW y»>l- "0 

^ ^ Lo 1 -1 -ij 


■ 2 with Mil- 


end transposed left latent vectors. 




m • 

3 -2 -1 -1 
_-5 5 0 5 ^ 


The right eigenvectors con be constructed froa (2.12) end (2.16) with be- 
ing the right eigenvector nstrlx 
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„Cl) 


X 

0 

1 


3 

X 

6 




.( 2 ) 


X,y<« 


-X 

-X 

X 

-X 


-2 

-X 

-5 

-3 


yO) 

^2 


y<»*X,y«> y<«tX,y«> 


. ,_(X) ^<1) „<2) <3>, 

• >'c2 >'c2 ' 


Th« trMUpoBed l«ft •iganvactora can ba conatructad froa (2.X3) and (2.18) 
with being the left eigenvector natrlx 


«L- 


(iXl+Aj) 

.( 1 ) 


T^-a) 


.( 3 ) 




,< 2 ) 

*2 


.(X) 


• 8 
X6 
3 
. 5 


X 

3 

-X 


X 

-7 

-X 


‘*cX *c2 *c2 *c2 ^ 


It should be noted thet the left eigenvectors end Xetent vectors ere 
reversed In order with respect to the right eigenvectors end Xetent vectors. 
The reeson for this reverse ordering will be discussed Xetera 

The coispenlon nuitrlx cen be generated froa the right eigenvector 
oetrlx or the left eigenvector aetrlx 


*c - «R •» "r" 


-T T 

V 


(2 a 21) 


where 
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J ■ 


10 0 
0 2 1 
0 0 2 
0 0 0 


0 

0 

1 

2 


In th« next aectlon It will be assuaed for simplicity that the (q^^-l). 
primary latent vectors do not have any generalized latent vectors and the 
q^th latent vector has gsnerallsed latent vectors although a more 

general treatment will require the structure of as given In (2.19). 


3. ElgetiproJ ectors and Latent Projectors 

It was shown In the previous section that the latent vectors of a 

regular lasdida matrix with leading coefficient - I are given as subvec- 

tors of the eigenvectors of the companion matrix constructed from A(X) . 

The structures of A^ and A(l) as given In (1.3) and (2.3) respectively will 

be assumed In this section. The computation of elgenprojectors and latent 

projections from the eigenvectors and the latent vectoxsvlll be given In 

this section and their relationship to partial fraction expansions of 
“1 —1 

[A.(A)J as well as [A(A)J will be given. The chain rule for the latent 

C 

projectors will be derived and an example will be presented. 

Zadeh and Desoar [1] have given the general partial fractlbn expansion 

of (XI-A^)"^ as 
c 


tA - (XI-A )*^ . I (3.X) 

' ' 1-1 1-0 (X-X^)**-^ 


where ^cC ™»m^ Is the multiplicity of the eigenvalue and q^ la the 
number of primary eigenvectors. The P^ ^ matrices will be called the elgen- 
projectors since they can be constructed from the eigenvectors of A^. The 
elgenprojectors are also the matrix residues which can be computed from the 
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uauAl foraulA 




. & 1 * / 

■ V 


1 d 


(k) 


'kl 


(x-v 


i.-q.+l 


(3.2) 


0.1, 




or frOB tho propor •oloction of tho right and loft oigonvoetora aa will ba 
soon latqr. 

Tho olgenprojoctora (or oatrlx roolduoa) oatlofy tho proportion of tho 
rooolution of tho idontity and opoctral docoapooition 


r 


I. 

^1.0 " ^anxan 


II. 

^i,0 ^1.0 ■ ^i,0 


III. 

^i,0 ^j,0 " °n«xnn 

1 f J 

IV. 

^.1+1 - ^.1 <V^i^> 

1 ■ 0*1.2. ... .a^^** q^-2 

V, 

^i,a^-q^^^c"^l^^ " °anXon 


VI. 

Ac “ I 0 ^ ^ 

i-1 ^ i-1 

'i.i 


C3.3) 


Mb 

Property VI doscribos tho opoctral doconpooition of A^eR and tho oot 
of oigonvaluoo ia callod tho spoctrua of A^. An additional proporty can 
bo addod for dofining functions of a aatrix which is obtalnod diractly froa 
VI with 




A-1 


1 d 


W 


\l 11 ^ 5 ^ 


H\) 


VII. 


(3.4) 
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fropartlti i~V «nd VII «r« glv«n in Z«d«h taxd Diioar [1] ^ 


3.1 Daflnltlon . L«t Wj^, Wj^eC*^"" ba formad fros tha «at p£ nonMllzad 
right and laftaiganvactora, and raapaetlyaly with - I. 

Lat «nd **• ^*** >f®®tangular watrlcaa (or vactora whan ■ 1) 

for with 





(q--!) <qJ 

• • • ^1 ^cl 


<-i>, 

^cl 


(3.5) 


whara It la aaaunad that tha flrat (^vactora In (3.5) ara tha prlnary algan** 
vactora and tha laft n^-q^ art ganaralliad alianyactora aaaoelatad to X^. 
Slallarly^ aaauna that 


'cl 


- [1 


(1) .(2) 


cl 


“cl 


(qi-1) 

'cl 


(■l) 

cl 


(q-) 

•cl 1 


(3.6) 


(1) (2) ^‘^1^ 

whara »***»*g;l^ 1“ (3*6) ara prlaary aiganvactora and tha vactora 

(o +1) 

- ^1 , the generallzad aiganvactora for X.. Tha ordarlng of tha 

cl ^ 

laft aiganvactora in (3.6) nuat ba aa ohown bacauaa tha onaa In tha Jordan 


blocka for X^ ara located on the aubdlagonal of due to tha tranapoaa 

operation In A . 

c 

Tha natrlcaa and ara tha right and laft alganvactor natrlcai (or 
nodal natrlcaa) with. 


*c - “r ■ “r O.n 

-IT 

where whim tha aiganvactora have bean properly nomallzad. It fol- 

lowa from (3.7) that 
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(XI-A^) - Wj^(Xl-J)W^^ - Wj^(M-J)wJ (3.8) 

thiw (Xl-A^)"^ “ WjjCXl-J)"^/^^, S«« thf uppendix for th« ntructuro of J. 

3,2 Thoorett . Lot b« « priBory •itonprojoetor of (3.2) for a 

distinct, nonrspsstod slgonvalus Xj^ of A^. Ths prlasry sl|«iiproJsctors srs 
glvon by 



T 

-VI 


^cl *cl 


(3.9) 


whars and arc noraalltad right and loft olgonvactora roapoctlvaly 

of A^. 
c 

Proof . Lot (Xl-A^)*^ bo ao glvon In (3.8) with P^^ ^ doflnod an In (3.2) 
than 


^.0- 


11a 

x-*-x. 


(X-X^)(XI-A^) 


-1 


- Y.z! - yfP 


‘1“1 


'ci *cl 


ao all diagonal oloaonto of (X-X^) (Xl-J)~^ will bo xoro In tha Halt oxcopt 
for tha 1th diagonal oloaonc which will be one. 

3»3 Thaoraa . Let P^ ***• priaary olganproj actor of (3.2) for 
a ropoated olgenvalue of A^ with - a^. Tha priaary alganprojactora for 
the rapoatad alganvaluas are given by 


P. « - 


1.0 


J-1 


cl ‘cl 


(3.10) 
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wh«r« ind ar« th« normllsad •igtnyttceori for 

froof. Tho proof followa frcm (3rS)^ Tho diagonal alaaaaca of 
(X-X^) (Xl-J)**^ will hava onaa only on tho dlagonala of tha Jordan block 
aaaoclatod vlth X^ in tha Unit and aaro alaavhara. Thla naana that aaeh 
right oiganvaetor paira off with a loft aiganvactor giving (3.10)# 

3«3.1 Corollary . Givan that y^^^ and a^j^^ ara not normliaad, than tha 
prinary aiganprojactora ara givan by 

(3.U) 

I j«i nj”' 

‘ci ^ci 


p . y'i ^ci «ci fi ycr«ei 
j-i J«1 


whara 

Proof. Aaauaa that tha nc^maliaad right and laft aiganvoctora ara acalad 
by and 3j raapactivaly, Tho nunorator of (3.11) than baeonaa 


- ojejycT'S’ 


wharaaa tha innar product In tho danoninator tamo ara 

\ *ci ^ci ®J®j*ci ^ci 
T 

ainci y^^^ - 1. Tho product in tha maarator and danoninator than 
cancali giving (3,10). It ia inportant to notica that a^^^ y^^^ - 1 bacaUaa 
tha prinary aiganprojectora aatiafy Proparty II in (3.3). 

3.4 Thaoran. Whan tha aiganvalua X^ ia rapaatad with nultiplicity n^ and 
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^ ^1 «l|tnv«ctori und g«ntrallit<l 

•ig«nvaetori. Th« priMry algtiiprojaefcor for la glvan by 





*ci 


(3.12) 


vhara tha flrat auBBatlon la ovar tha aat of q^-1 noraallxad prlaary algaii- 
vactora md tha aacond auaaatlon la ovar tha qth prlaiary alganvactor and 
tha aat of noraallsad ganarallzad alganvactora. 

Proof. Aa baforat can ba wrlttah aa In (3.8) with J. Utilizing 

mmmimmtmtamm Q 

tha natorial mm giyan In App«ndlx A, the prlaary aise^pcoj actor la 


^.0 ■ (>■ «) 


which aftar dlffarantlatlng (■^-c^^)tlMa and taking tha Halt glvaa 


1,0 


'‘Ai. 


(3.14) 


and aubatltutlng for and aa daflnad In (3.5) and (3.6) gives (3.12). 


3.4.1 Corollary . Lot y^^^ and z^^^ be unnoraallzed right and left eigen’* 
vectors for tha algcnvaluo which has nultlpllclty and < a^. The 
prlaary alganproj actor for X^ Is given by 




■•■A 


f ypl ^cl . rl ^ci *ci 

Jj ^ 


3-1 


3^,1 



(3.15) 
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iriun y^i^ 


Froof i Th« proof of this corollary followa fro» tha proof in Corollary 3*3 f1 


3*4*2 Corollary . Tha nonulliatlon coafflclanta ri^* for Cha auMutlon 
ovar tha qth prlaary alganiMetor and tha ganaralliad aiganvaetora la glvan 
by 


** • * * * ^ 


‘* 1»*1 


(3*16) 


Proof * Tha proof of thla corollary follova liMdlatoly fro* tha daftnltlon 
of ganarallsad right algaavactora aa givan in (2*18) and froa tha guaral- 
iiad laft alganvactoxaaa follova 


- (*1+1) 
‘cl 

• ‘cl 

m (q>2) 

(*1+1) 

*cl 

-‘cl 

■ <v 

T ^ 
*cl 

■ *cl 


and taking appropriata Innar product tha follovlng raaiilt la obtalnad 

T 


’^1 *ci ^cl *cl ^1*' ^cS 


for <1^ 5. 3 1. “i — b ^ auch that j|+k ■ than (3.16) la provan* 


3.5 Thaoran * Lat y^|^^ and ba tha right and laft alganvaetora for 


tha aiganvalua 1. of aultiplicity a. and q. < ■.« than thara will ba a.-q^ 
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*1 


g«iMiralli«d alganprojactora •■■ociaud to and givan hy 

y_ .%T 


IfA 




0.17) 


whara I ■ lp2i • • • 

Proof . Thia thaoraa la provan by conaldarlng tha atructura of 

aa givan In Appandlx A» and by carrying out tba algabrale oparatlona of 

(3.2). 


3.5.1 Corollary . GlVan that tha aacondary alganprojactora ara to ba con- 

atructad froai tha aat of unnoraalliad alganvactorai tha ragulrad noraallsa- 

■i.qi 

tlon factor la » l.a. 








*cl 


A*l*2p » . . (3.18) 


4 V 

Proof . The normalization coafflclanta ara conatanta ovar tha qth 
primary and ganarallzad alganvactora of Xj^, tharafora^aach vactor product 
In (3.17) will be nomallzad by the aane conatant. 

The algebraic theory of elgenprojectora la uaeful In developing the 
algebraic theory of latent projactora (or matrix raalduaa) of a lambda 
matrix. The partial fraeft^^on expanalon of tA(X)]’’^ la 


[Awr^- I 0.19) 

lil t-O 

A A 

where P^ la a latent projector. It la known that P^ can be computed In 
tha uaual manner with 
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ila 

X-^A, 




1 d 


<k) 


(A-A^) 


^ ^ iA<A)r^} 


(3.20) 


for k ■ 

Tho dovalopoant on •Igonvoctoca, olgonprojactora and latont voctora 
will now ba uaad to davalop tha algabraic thoory of latant projactora aa 
dafinad In (3.19) and (3.20). 

t 

3.6 Thaoraa . Tha laabda natrlx [A(A) will ba glvan In tha uppar right 

a^B block of tha invaraa of A^(A) an ahown 

c 

lA^(A)]"^ - 


IA(A)] 


*•1 


(3.21) 


providad that A(A) la ragular and Aq ■> I. 


Proof . Tha proof of thla thaoraa followa dlractly from tha Invaraa of 

A^(A) whara 
c 


A^(A) 


Al -I 

0 AI 



0 . . . 
-I 

. . . . 

. . . . 



0 

0 



/N 

3.7 Thaoraa . Lat P^ q danota a prlaary latant proj actor of A(A) for tha 
latant root A^ of aultlpllclty a^ and q^ ■ a^. Tha prlaary latant projactor 
la glvan by 
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ft . yi 

^,0 *1 


(3.22) 


M) -<J) 


whar* and ara nonullzad rlgltt and laft latant vactora raapactiva- 

ly of A(X). 

Proof . Tha natrlx A^(X) can ba daflnad by 


[Ac(»r^ - wj 


and front tha pravloua thaoraa 


[A^(A)1 



1 

r< 

< 

a 

a 

a 

[ 


Y 

1 _ 

• f • a 

4 • 9 9 

■ 

s 

s 


ass 9 


a 


(M-J)"^[....z’^l 


with Y " [Vj^,y2» ” * ^ ^r^* latant vactora in 

and ara In aaquantlal ordar since thara are no generalized latant vactora. 
Using tha definition c/f an ai|k)anproj actor from (3.20) and 


[A(A)r^ - Y[Xl-Jl”^z'^ 


gives the desired results when the individual vactora of Y and Z^ ara con- 
sidarad. 

3.7.1 Corollary . Tha laobda matrix [A(X)]'*^ aa defined in (2.3) ia given 
by 


[A(X)1“^ - Y[XW1“^Z^ 


(3.23) 
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wh«r« Yec“^"“ «nd ZeC*^"“ «r« th« matrlc** o£ right and Xaft lataac vactora 
raapactlvely. Equation (3.23) will ba daalgnatad aa tha canonical fom for 
tha Invaraa of a lairisda aatrlx [3j. 


Proof s Tha raaulta atatad In thin corollary warn obtalnad la tha pravloua 
thaoraa. 


3.8 Thaoraa . Tha primary latant projactora of AtD for tha latant 
root of Bultlpllclty and q^ < ara glvan by 



whara yp^ and ara tha noraallzad right and laft latant vaetora ra- 
apactlvaly. 

Proof . Tha proof of thla followa froa a datallad analyala of (A(X)]~^ aa 
daflnad In Corollary 3.7.1. 

3.9 Theorem. Let ba a primary algenproj actor aaaoclatad with 

XpU 

an elgenvalua of multiplicity with ■ a^. Tha a><B block aatrlcaa 
In tha laat nr'coluana of g ara glvan by ^ " 0,1,2, ... ,n-l 

whara P^ g la a primary latant projector. 


Proof. Whan tha Invaraa of A_(X) la tak«i tha a laat coluana will ba glvan 

c 

aa followa 


(A^(X)f^ - 


lA(X)f^ 

. X[A<X)f^ 

. X^[A(X)1"^ 

a a 

a a 

a a 

. X"“^(A(X)1“^ 


(3.25) 
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but th« •lg«BproJ«ceor ^for is dtflnsd 


' i.o ■ 


(3.26) 


than applying tha last aquation Into (3.25) glvaa 


1,0 


• a 


a a a 


a a 


a a 


a a 


'1.0 

Vl.O 


^1 ^1,0 


(3.27) 


which proves the thaoreD. 


.mnXmn 


3.10 Thaorem . Let P^ ^eC be an eigenprojector for of aultlpllclty 
with q^ primary elganvectora and n^-qj^ ganarallzed elgenvectora. Tha 
mxa block matrices In tha last m columns of ^ are given by 


1.0 


a a a 


a a 


• a a 


a a 

a a 

a a 


a a a 


'1,0 

^l^l,0'*'^Vl,l'’‘^1.2 


j,n~li . in“l» . n—2p ,^n— lv»n—3p , 

Xi ' 1 ^'1 ^i,r ^ 2 ^^1 ^1,2^* • • 


(3.28) 


Proof . This theorem Is proven by considering the eigenprojector P. ^ de- 

l,u 


fined as 
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<«l-qi) 
. 1 d * ^ 

^1.0 " X-+-X^ ^ I 

dX 


a-\) ^ ^ EA^(x)i "•) 


(3.29) 


when this property la applied to (3.25) » then the mKa block .aatrlcaf In the 

laat coluana of P. atructure aa given In (3.28). 

i,u 


3.11 Theoren . Let be the latent projactora eaaocletad to the 

regular lanbda matrix A(X) with A^ - 1, the aet of latent projactora aetla- 
flaa the following propertlea 


X. 


I - 

I '^.0 


- 0 


bXib 


r i 


“■ ill io *" ^ ^ ”■* -"‘■'i 

«!• T ‘iv»i »-w 


If ell then all the generellxed latent projactora ere aero and the 

above propertlea alopllfy to 



j “ 0,... ,n-2 


(3.31) 


Proof. Let P. £0““*^““ be the primary elgenprojectora of the companion 

r 

form A aa given in (3.28). when the reaolutlon of the Identity ^ P. _ ■ I 
^ !■! * 
la applied, the mxm block matrlcea in the laat column aatlafy the propertlea 

given In (3.30). If ell - m^ then P^ ^ have a almpllfled atructure ea 

given In (3.27) then propertlea In (3.31) can be eeally obtained. 
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Proparty 1 In (3.30) will ba callad tha raaolution of tha zaro owitrlx. 
Tha abova propartlas have baan darlvad whan tha lanbda natrlx haa tha Idan- 
tity as tha landing matrix coafflclant, but for a genaral ragular lambda 
matrix tha abova propartlaa alao hold whara Fropartias 111 and V ara nodl- 
fiad having In thalr right hand slda. 

It haa baan aaaumad; In tha thaorama on latant projectors that tha 

m 

latent vactora ware nonaallzad with normalization ■ 1. The abova nor- 
malization la based on complete knowledge of tha algenvactora which may not 

ba available. It la tharafora necessary to compute tha normalization 

1 k 

factors from the latant vectors. Suppose that la as given aarj,lar, 
than for normalized algenvactora 


^1 *ci ^ci * 


(3.32) 


This normalization factor muat ba applied to each lumormallzad latant vector 
In tha latant projector calculation, 

Tha primary alganvactors of ara given In (2.12) and (2.13) tharafora, 
tha normalization factor, la 


j.j 




Ji) 

^1 

X y^^ 

'^ 1^1 


X“-Xy<4) 

Ai 


.(3)^ (3) 

*1 dX 


for 3 < 


(3.33) 
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Th« nonMll;i«tlon f«ctora for th« latrat vactora in th« chain rule ar« 
Bora coaplax chan tha nomalizaClon faeto)ca for tha prinar^ latanC vaetora. 
A laana will ba givan for that caaa. 

■l^l 

3,12 Lfa . Lot J^eC ba tha Jordan block aaaociaCad with tha laeant 
root and lot and a^^ ba tha laCant vaetora. Tha noraaliaation 
factor la given by 


J.k 


!. 

t+a4p-j+k+g^ 


.<t) 


JL 

at 


.(•) 


dX 


A(X^) 

(JT” 


-<P) 

^1 


(3.34) 


whara X^ haa ganaralizad latent vaetora and irisara <1^ £ t £ j and 

•ll 1 P 1 


Proof . Tha proof of thla luana followa froa Corollary 3.3.1 and tha corra- 

apondlng ralatlona for the aigenvac^ora and latant vaetora. Aa •t'ttad 

■l.5l ■4-i»q4+i 

Corollary 3.4.2, tha nonaallzatlon factora ... ■ 

‘*l*“i 

aro oqual. 

Tha gonerallzod latant projactora aatlafy a chain rula which would bo 
axpactad olnco thoy arc conatructod fron tha ganarallza latent vaetora that 
aatiefy a chain rule. 


Thaoran 3.13 . Lot ganaralizad latent projactora aaaocl- 

atad with tha latant root X^. Tha chain nila for tha ganaralizad latant 
projactora la 


■i‘K^ i -0 


J " 0,1, . . . (3.3S) 


where A(X)i:R**^* la a regular laabda aatrlx. 
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m^Bn 

Proof * Lot P^^jfC bo tho olgonprojoctori of tho coapanlon fora A^. 

Uolng oo doflnod In (3.2) It la oaoy to show that tho hKb block ratrlcoa 
of tho laat ■ eoluana havo tha follovlng atructuro 


l.A 




>h“lp 


(3.36) 


whan tho structure of ^ la used with Proportlos IV and V In (3.8) then 
the results in this thoorom are obtained. 

Tho results In the last theorea also hold for regular lanbda matrix 
where matrix leading coefficient Is different than the Identity matrix. 

Tho following example will be given to Illustrate some of tha computa- 
tional aspects of the material. Lot A(X) be defined as 


A(^) 


1 0 
0 1 


+ 


•4.0 2.0 

-0.2 -3.0 


X + 


3.0 -1.0 

0.2 2.6 


with latent roots Xj^ ■* 1 of multiplicity 1, and X 2 " 2 of multiplicity 3. 
The latent vectors for the lambda matrix are given by 


and 


y<« y‘«l - 




1 3 -1 -2 
0 1 - 1-1 

3 -2 -1 -1 
■5 5 0 5 


4 ! 




II 


ll 


% b 
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Sines nullity A(2) I, thsra is only ons priasry Istsnt projsctor for 


^2 " ^ * 


Ths Istsnt projsctor for * 1 esn bs eonputsdi dirsctly fron ths Istsnt 
vsetors y^^^ snd and dA<l)/dX or 

.T 


„(1) -<1) 

a *1 


1,0 


1 

5 


-3 5 
0 0 


gy dA(i) (1) 

*1 dX ^1 

Ths sscond priaary latsnt projsctor ^2^0 by 


(2)* (3) (D* 

^2 *2 ^2 *2 ^2 *2 _ 1 


‘ 2.0 


.3.1 


3 -5 
0 0 


3 1 

vhsrs n,* is coaputsd from (3.34) with 


3.1 . < 3 )* dA( 2 i (1) + 1. .(2)* A(2> (1) X_ ,<X)* Am. (1) . .5 

^2 " *2 ^X^ ^2 ^ 21 ■ jj^2 ^2 ^ 31 * ^2 ^ 


It can bs vsrifisd that “ hJ*^* Ths |snsralissd latsnt projactors 


ars 


„(1).C2)^ <2) (1)^ 

^2 *2 ^2 *2 _ 1 


2.1 


,3.1 


2 5 
0 5 


and 


2,2 „3.1 


(1).(1)‘ 

2 *2 . 1 

I 


3 -15 
1-5 


Ths aigsnproj actors for ars 
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and 


1.0 


2.0 


2.1 


2.2 



1.6 

- 3.2 

- 0.6 

1 


0 

0 

0 

0 


1.6 

- 3.2 

1 

o 

m 

1 


0 

0 

0 

0 





— 







- 0.6 

3.2 

0.6 

-1 

B 

0 

1 

0 

0 


- 1.6 

3.2 

1.6 

-1 


0 

0 

0 

1 


mm 










- 0,4 

- 3.2 

0.4 

1 


0 

- 2 

0 

1 


- 1.4 

- 2.2 

1.4 

-1 


- 0.2 

- 2.6 

0.2 

1 





— 


...» 





- 0.6 

4.2 

0.6 

-3 


- 0.2 

1.4 

0.2 

-1 


- 1.2 

8.4 

1.2 

-6 


-0.4 

2 . a 

0.4 

-2 


n 

u 


I } 


iriiara it !■ notad that tha 2«2 uppar right blocka ara tha latant projaetora. 
Tha aigenproj actors wara eoaputad by using ah ^ w“^ with 



as eonstructad from (3.5) and tha shifting natricas ara 


( 

i 
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■h r 


1|0 


•h r 


2,0 


10 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 


0 0 0 0 
0 10 0 
0 0 10 
0 0 0 1 


•h F, , - 

*»* 


0 0 0 0 
0 0 10 
0 0 0 1 
0 0 0 0 


«h r 


2.2 


0 0 
0 0 
0 0 


0 

0 

0 


0 0 0 


0 

1 

0 

0 


Th« structurs of the thlfting matrlcea h«a alao appeared In [10], 


The partial fraction expanalon for [A(X)] la 


1 1 

-3 5 

1 

3 -5 

t ^ X 

2 5 

" six-iy 

0 0 


0 0 

5(X-2)* 

0 5 


5(X-2)' 


The chain rule for the latent projactora ••tiafy tha aquatlona 


I. A(2)P2^2 ■ ° 

A(2)P + P A " 0 


II. 
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given In (3.35)* 

The ranulta given In thin chapter hold only for luibdn lutrleoo that 
nro rogulari l.a.,AQ la nonalngular and a full aat of latont vaetori axiata. 

4. Conclualon 

Tha theory of latent projactora have bean praaantad and their uae In 
the conputatlon of the Invaraa of A(X) and Ita railduaa have bean deacrlbad. 
Tha latent projactora have apaclflc propartlaa and are ralatad to tha algan- 
projactsrfl of the conpanlon fon A. obtained from A(X) . The concept of 
gimarallzad latent projactora haa bean developed which are conatructad fron 
tha generalized latent vectors of A(X). Tha chain rule for tha generalized 
latent projactora haa bean developed and deacrlbad. The work on eigen- 
pro jectora, latent projactora and laabda aatrlcea preaented In thla chapter la 
not exhauatlve aa mny other queatlona nuat be reaolved; for exanple, the 
caae of alngular A^ haa not been addreaaed* The next chapter raaolvea the 
problem of alngular A^ • 


CHAPTER 3 




THE EIGENVALUE-EIGENVECTOR STRUCTURE OP OBSERVER STATE SPACE PORKS 
AND THEIR RELATIONSHIP TO GENERAL LANBDA KATRICES 

Thlf ch«pt«r itudlti th« Algcblriie theory of latent rooti end latent 
vectors of e row-reduced fora of a laabde matrix A^ X ) ■ AqX"4‘Aj^X^~^ 4’... 
,.i-fA^« In the preyloua chap tare Ap ^ waa aaauaed to be rigular» i.e^i 
dat Aq ft Oj however, thla aesunptlon Is very restrictive particularly When 
A( X ) arises as the left Inverse matrix In a left polynomial matrix descrip- 
tion for a dynamical system. In this chapter the regularity assumption Is 
replaced by a weaker assumption, namely that Aj^.CX) Is now reduced. The 
latter asaumptlon Is shown to be without loss of generality In the context 
of dynamical system modeling. The development utllltlzes the duality 
between lambda matrices and observer state space forms to reveal the algebraic 
structure of the latent vectors of row-reduced lambda matrices. This, In 
turn, exposes the algebraic structure of the eigenvectors of observer state 
forms. 

This chapter along with the material presented in Chapter 3 gives the 
general algebraic theory of lambda matrices and their relationship to state 
apace forms. The mathematical development In the chapter supplements the 
material In the previous chapter. The co^>utatiOnal procedure for the latent 
vectors, eigenvectors, latent projectors and elgenproj actors for the singular 
case Is given. 
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Hultlv«r>i4bl« •ytUMi, luch as flnlta-alaMinl; nbdala of ylbratlni oya- 
toM, eoatvol ayatoM^ and laxfa acala ayatmap can ba fonailatad in taraa 
of aaeond-and hi|hai^*ovdar iutrlac dlffarantlaJ. aqua£lona of tha foni 


A -^ + A, 
0 dt" ^ 


dt 


n-i 


+ 

II 



<¥ in 

U 

a.i) 


whoxa x(t)cR*p and u<t)eR^« Such a aodal la of tan eallad 

a loft BACXlx dcacrlptlon for tha dynaalcal ayataa. 

If tha Laplaca txanafora la tafcan of Eq. (1.1) i tha follewlsg aatrlx 
polynoaial appaara 


ACa) ■ Aga“+Aj^a®"^ + + A^ 


( 1 . 2 ) 


on tha laft alda of (1.1). If a la raplacad by Xp than (1.2) iaflnaa a 
lanbda aatrlx [1], alao eallad a natrlx polynomial. 

Our objactlva hara la to Invaotlgata tha latant roota and tha latant 
vactora of tha polynomial natrlcaa. Ear liar work on thla topic haa baan 
daacrlbad by Lancaata:!? [ 1 ]. Lancaatar and Wabar [4]p Dannlap Traub and Wabar 
[S]p Lancaatar (dj. and Danmanp Layva~Ramoa and Jaon [U]. In all of thla 
work it haa baan aaauaad that A(X) la raaular . l.a. , dat A^ ^ 0. Thla 
aaaunptlon parmlta tha modal in (1.1) to ba axpraaaad in block companion 
atata apaca form. Tha aaaunptlon that A(X) la ragular la too raattictiva 
In the contaxt of modala of tha form In (1.1). In general, regularity can 
ba achlavad by Introducing additional roota, l.a.p at tha origin. Bowavar, 
thla nagataa ona of tha prime motlvatlona for tha work In [3]-[5]p namely 


that of achlovlng nunarlcally officiant ptocaduraa for handling nigeinvalua- 
oiganvactor probloao ariaing froa ayatoaa of Kha fom givan in (1.1). Tha 
ganoral caao will bo traatad in thia chaptar without tha naad. for introducing 
additional rooto. 




OWGlNlJt »« » 

0^ fOOR QUALfTT 

2. Row-redueid Lambda Matrlcaa 

To onauro oxlitonco and unlquanoBB of tho aolution to (1.1) » It la ra-> 
qulrad that dat A(X) ^ 0. Tha following daflnltiona ara intreduead, 

[131 

Daflnltlon 2.1 A nonalngular polynonlal matrix Aj^d) !■ ttld to ba tcM- 
raducad If tha matrix fomad from tha coafflclonta of tha hlghaat poWara of 
X In aach row la nonalngular. 

It followa that any row-raducad polynomial matrix Aj^(X) can ba wrlttan 
In tha form 

yX> - S(X)Dq + L(X) (2.1) 


whara 


s(X) 

- dlag[X ^,X ^ X ®] 

(2.2) 

L(X) 

- fCX)*^ L 

(2.3) 


V--1 V -2 I 

X ^ X ^ 1! . 


T 

‘Tv^ V2 

:x ^ X ^ ... 

1 

« 

1 ! 

a a a A j 

'i'(X)^ - 

• 

a 

• ^ f k .k ft 



|X* X® ... 



(2.A) 

and Dq, L 

are matricea of raala with 

, 

dat 

“o “ 

(2.5) 


A kay fact from tha thaory of polynomial matrlcaa la tha following: 
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L«— 2.1 Any nonalnpilnr laabd* Mtrix A(X) can ba tranafonad to roit- 
raduead fom Aj^(X) by laft wiltlplleation by a unlaodular aatrlx U(X). 

Proof: Saa [li P* 271 or tal. ^ 

Nota that operating on £q. (1«1) by U(X) doaa not change the number of 
latent roota alnce det U(X) la a nonalngular conatant matrix Indapendent of 
X. Thua,'wa aee that A(X) can be tranaformed to the row-reduced form 
without loaa of generality . However , It la vary raatrlctlve to aaaume that 
Aq In (1.1) la nonalngular alnce this la tantamount to aaaunlng that the 
row Indices to are all equal. Thit latter situation la unlikely to 
occur wilaia additional roota are artificially introduced. Of course, the 
case det yi 0 is simply a special case of roi^reducodnoaa in which 
'>1 ■ 'V' 1, Je[l,2,. 

m 

Theorem 2.1 . The number of latent roots of A(X) la precisely ][ V. where 

1-1 ^ 

la the row-degree In the row-reduced form of Aj^(X)> 

Proof: The proof follows directly from (2.1) as dat ^ 0. WV 

By standard argumenta, [L3], [£<4, any row-reduced set of differential 

aquations, as given In (2.1)^ can be associated with a state space model 

» 

In the observer canonical form, aometlmma called the flat canonical form 
pLf] . The system matrix has the following form: 
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1 0 M . 0 I 
. • • • 
0 . • • 
V’ , * I 
* ’ .0 • 

** ll 


'pI 


F ■ 




0 M 


1 0 

. ' 

0 * 


0 •* 


I 

I 

I 

I I 


( 
» 

I 

I 
» 

I 

I 
I 

I 

I 

I 
I 
I 


1 0 0 

* » . 

0 • 

’’ • -A 

• , 0 

• « 

• 1 

0 


( 2 . 6 ) 


wharc [LDg^] ^ danotas tha 1th colu» of LD~^ and Dg and L ara aa In (2.1) 
and (2.3). 

The above companion form la tha appropriate generalization of tha block 
companion form uaad alaewhara, (aaa [a] to [5]). Nota that tha block com- 
panion form requlrea that dot Ag 0, wharaaa the form In (2.6) neada only 
tha weaker a8ausq>tlon that dat Dg 0 and there la no loaa of generality In 
tha latter aaaumptlon. In tha aaqual, * 1 will be taken alnca If thla 

la not tha caaa, than one can almply multiply by Dg^ throughout. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


43 


3 . Eigenvectors and Generallzad ElgenvecnogB 

In this icctlon, the definl^inni for elgenvnluesi right end left eigen- 
vectors, end the cheln rules for the generellxed eigenvectors of e reel 
■etrlx F will be given. 

The following definitions ere stsnderd for e reel witrlx F, (see [16]i 


Definition 3.1 A sceler is en eigenvelue of F if 


det[X^l-F] - 0 . (3.1) 

The eigenvelue X^ hes oultlpllclty n^ if 

Q n. 

detlXI -FJ - n (X-X.) (3.?,) 

i-1 ^ 

where a is the nunbsr of distinct eigenvelues. 

Definition 3.2 A vector is seld to be sn eigenvector of F eesocleted 
with X^ iff 

(F-X^I)y^^^ - 0 (3. 3) 


(2) (3) 

end e set of vectors {y^ ',y^ 
eigenvectors essocietsd to if 


(h.) 

,y^ } ere seld to be the generellxed 


(F-X.dM*"^ -0 k - 2.3....,h, (3.4) 

^ (h >1) ^ 

where h^ is such thet y^ 0 is the only solution in (3.4) . 
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Tha gcnarallzed aiganvactora , aa givan in (3.4), can ba alan obtialnad In 
tha following fona 

* « 

(F-l^I)yJ*‘^ - k»2.3,...,h^ (3,5) 

k 

and tha ganaralizad aiganvactora aaaociatad to y^^^ font a chain whara h^ 

ia the length of the chain. It ia posaiblo that tha aana eigenvalue produce 

aata of differenc chalna aaaociatad with other aiganvactora but tha 

number of different chaina aaaociatad to the eigenvalue la equal to tha 

nullity of [F»X.I]. Tha length of tha longaat chain la called the index of 

. g maxh. 

X, and is such that \|>(X) ■ H (X-X. ) ia the minimal polynomial of Fi 

1-1 ' 

Tha expression In (3.5) is called the Jordan chain of generalized right 
eigenvecto rs. 

From the eight eigenvectors the modal matrix can be constructed by the 
collection of the eigenvectors and their chains where the chained aiganvac- 
tora appear in a sequential taanner aa Wj^ ■ ,. . . iy 2 ^\y 2 ^\ . . . » 

yi^\yl^\ • . •] • and Wg is called the right eigenvector matrix. 

A similar definition of tha left eigenvectors of F can be given where 

(F-X^l) is replaced by its transpose in definition 3.2. Tha sat of vectors, 
(2) (3) (h,) 

(z£ '^,...,z^ } form a chain of left ganerallzed aiganvactora asso- 

ciated with tha primary eigenvector z^^^ . Tha chain can be obtained from * 

(F-Xj^I)'^zJ*'^ - k «2,,,,,h^ (3.6) 

where h^ is the length of tha Chain and z^^^ is the primary left eigenvector. 
Also, tha collection of all the left aiganvactora forma the left eigen- 
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Mtrlx vh«r« tha chalnad aiganvaetora appaar In a aaquantlal naanar 

fyd) _(2) (1) „(2) <1) (2) , 
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4, Latent Roof and Ljtf nt Vactora o£ hMwMm Matrices. 


Anelogoue definitions to those In Section 3 for s roie-reducsd Isnbds 
■etrlx Ap^(X) will be given. It will be shown that there Is an Interesting 
duality betwasn the definitions of Section 3 and those for the correspond- 
ing lanbda matrix A^(^) m given In (2.6). Without loos of general- 

Ityi It Is assunad that Dq ■ I. To simplify tho notation in the proofs, It, 
has boon assuPMd that 9^2. The extension to a>2 Is along the same lines. 

Definition 4.1 Let Aj^(X) be es defined In (2.1) and let X^eC be a scalar, 

then Is a latent root of If det * 0* Now hae multiplicity 

n. If det A(X) ■ II (X-X.) where Q la the ntmber of distinct roots. 

^ 1-1 ^ 

Connection 4.1 The determinant of Aj^(X) Is equal to the determinant of 
(XI-F), where F Is as defined In (2.6) with 


det Aj^(X) - det(lX-F) 

^ “d 

- n (x-x.) ^ 

1-1 ^ 


(4.1) 


Proof: Let Ap^(^) be es defined In (2.1) and represented by 


Aj^(X) - 


11.. . 12. vx 

p (X) p (X) 

21.. . 22... 
p (X) p (X) 


x'^^+ej^x'^^ -f + 

e21^V2-l 4 . • + 

• 1 V- 


12 

a“x ^ + 

.^2 22,^2“^ . 

X "^ej^ X + 


. . . . "^ a 


. . . + a. 


12 

‘’l 

22 


(4.2) 


ess 
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L in (2.3) can ba aaally obtainad aa wall aa tha obaarvar canonical fon If 
than (Xl-F) can ba rapraaancad by 





If (Xl-F) la aultipliad by the following unlnodular aatrlx. 


> 


I ’ 

than 


1 0 ,. 

X 1 . 


U(X) ■ Block DlagV | ‘ . . 

I *vi ■•*.'.’.0 

^ . ; . . . X • I 
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U(X) (XI-F) 









•J 

* •, 

• 4 • 

0* • • •*( 

1 

a 



a 

• 

1 


a 

0. \ '0 

a 

a 

« 


a 


' 4 • 

4 

9 

• 

• 


p“<>) 

1 * • 

: '-.-i 

p“<» 

• 

a- 

a 



* % 

0 • • » • *0 


• 

0 « * « ii ( 

• • 9 • • a 








Pu 

Vi 

1 

)•••*• »g 
• 

* 

“1 0 *■ * *c 

* » 

'4 

1 

» 

1 



a 

a 

‘ \ 

1 

: 

a 


. 


0 * 0 

a 



a 

% • 


21/ix 

Pj (A) 


* 1 
• 

p“(X) 

• 4 * 

; 4-1 

a . 

P (a) 

n 

/ 

• « 

p»(X) 

0 * . . . t 


V , 


J 


G(X) . (4.5) 


Howavar, by avaluatlng |g(X)| by ainora on tha onaa, It la cl^iar that 


dot(Xl-F) ■ dat 


■p“a) 

p‘’<X) 

p“‘(X) 

P^^X) 


dat A(X) . 


(4.6) 

m 


Daflnltlon 4.2 Lac Aj^(X) hava latant roots X^ than tha prlaary right latant 
vactora danotad by y|^^ and tha prlaary lalt latant vaetora danotad by 
satisfy tha ralatlons 

and 

(4.7b) 


raspactivaly. 
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W« nw axplora th« ■tnictur* of tho •Igonvoctori of F uioclotod with on 

•Igonvaluo and ttion rolato thoia to tho lotont voctora uioelatod with • 

(li 

Itom (3.3) tho ttigonvactor iatloflos tho rolotion (T-X^X)yQ^ ■ 0 and 
by oparatlng on t.hli axpraaalon by -U(X)t ao daflnad In (4.4), tha foliow- 
Ing axpraaalon la obtalnadt 

® 

whart G(X) la aa daflnad In (4.5). Nota that tha polynoplala Pj^(X)» 
j - ara tha Tachlnhautfar polynonlala [16] aaaociatad with p^^(X) 

daflnad aa followa 




fi 1 

n 


r 


e • • 

• • » 

• 


1 

• 


2 • . 

X* . . • • 

• . . • » 

• 

» 

a 

f 


,11 

1 


■ 

a 

. t , * * 

a ♦ • 

f 

• 


f 

. p^‘(X) . 


‘ ■ 1 ■ 

^ * • • . 
V, 4^ v.-l • 

[ X ^ X ^ ... ’ X^ X 

• " 
0 

« 

1 


a 


and almllar axpraaalona can ba found for tho othar polynoalala In (4.5). 

Tho alganvactor can now bo axproaaad aa a function of tha latent 
vector y^^^ In tha following connoctlon. 

Connactlon 4.2 y^^^ bo aa dafinad In (4.8) and lot y^^^ ba tha eorre- 

apondlng right latent vector aaaoclated with the root X^^^, then the eigen- 
vector and latent vector are related by 

yo“ • 


<4.10) 
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vh«r« U(X^) !■ M d«fin«d in (4«4) and M la flvan by 






•i' 


,11 

\-X 


Vi 




.22 


.21 


U "2 


Vn-1 


22 

V«-1 


(4.U) 


Proof ; Tha axpraaslon glvan In (4 #10) la obtalnad by noting In (4.5) that 
tha firat and tha (v^-fl) alananta of tha alganvactor form tha latant 
vector y|^^ . Than by ualng al^>la aubatltutlona to obtain tha raat of tha 
alananta of the alganvactor aa a function of tha alananta of tha latent vec- 
tor, tha axpreaalon (4.10) la obtalnad where tha corraapondlng Tachlnhauaar 
polynonlala aa glvan lu (4.9) have bean uaed. 77V 

Alao tha loft eigenvector of tha obaarvar fom defined aa 

T (1) 

(X^l-F) Zq^' ■ 0 can be obtalnad aa a function of tha left latent vector 
z^^^ of tha row-reduced fom Aj^(Xj_) : 

Connection 4 .3 Let z^^^ be the left alganvactor of F and the corra- 
apondlng laft latent vector aaaociatad to the root X^, than the alganvactor 
and latent vector are related by 


original page m 

OF POOR QUAUTY 


-CL) . 




(4.12: 


wharf ¥(X) la as glyan in (2.4). 


Proof: Conaidar (Xl-F) aa glvan In (4.3), than tha laft aiganvactor la ob- 
talnad by (X^I-F)^Zq^^^ 0 and using tha axpraaalon for *qJ^» than 


11 -11 


'*t«r 

-1 X 0 • • • • 0 

0 * . • . •• : 


11 -21 


x"^ 

^1 


a“ 

•1 •* 


0 -1 • X. 


* .12 ...22 22 


X+*‘* m** 

^ 1*1 *2 • • • 


-1 X. 0.. . 0 

. 1 V • 


0 

• • a 

• • 


•. 0 


0 » . . Q -1 * X 




In tha abova axpraaalon, tha flrat and (Vj^+1) rows will contain tha laabda 
matrix A^^(X) aa daflnad in (4.2), tharafora ■ 0 aa axpactad. V7V 
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It it inttraitlng to notico In (4.12) that tht tloaonta of tha laft 
latant Victor appear on the and locatlona of the laft alganvactor. 

In tha aame aMuuar aa In paction (2), a aat of ganaralliad latent Vae- 
torp aaaoclatad vith the latant vector axlat and they aatiafy a chain 
rula, which Ip given in the next theorex. 


Theorea 4.1 tat A^C^) he defined aa in (2.1). A aet of right latent vactorn 

,(1) ^(2) ..^*'1^«BX1 

'l '^1 


7J~' iYj ' * fom a right Jordan chain aaaoclatod with the 


.( 1 ) 


latant root and the priawry latent vector y^ The chain ia given by 

ir 


2 

.a )v“> + ,<*■“ + i- * .<*-*> + 

W''i * dT^>'i +21 3 ’'i +••■ 


. . . + 


( 1 - 1)1 




Vl 


1 * 2 1 i . » |h 


(4.14) 


( 1 ) 


where h^ la the length of the Jordan chain. The vectora y^ for l<l;^hj 
are called the generalized right latent vectora of the latent vector y^ 


( 1 ) 


( 2 ) 

Proof ; The generalized right eigenvector yQ^' can be obtained f roai the 
equation (T-X^I)yQ^^ ■ y^j^^ which when Bultlpllad by -U(X^) aa defined in 
(4.4) givea the following expreaaion: 




(4.15) 


Uelng connection 4.2 and conalderlng the and hhe above 

expreaaion, the following reault la obtained 




( 4 . 16 ) 





- 0 . (4.17) 

(3) (3) 

Coasldarlng the generalized eigenvector as obtained froa (^'‘X^I)yQ^' ■ 
yQ^\ than a alallar axpraaaloh to that In C4.16) la obtained with 

- 0 . (4.16) 

yiP la needed In the above expreaalon and can be obtained from (4.15) where 
'01 

the eleiMinta of y^^^ are located In the flrat and (V^+l) aleoanta. We then 

( 2 ) 

uae back aubatltutlona to obtain the reat of the eleaenta of y^^' aa a fimetion 
of y^^^ where the following expraaalon la obtained for y^^^ 
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Th« Mcond aatrlx on fch« aocond tarm la not U(X) but la a ahlftad var- 
oion of U(X) and will ba callad ah U(l) %diora ah atanda for ahlftlns. Thla 
ahlftlng Mtrlx haa tha Intaraatlni proparty that 


1 # 1 ^^) 


(4.20) 


Mow aubatltutlng tha valuaa of y^^^ Into (4.19) and ualng tha ahlftlng 
proparty, tha axpraaalon In (4.18) can ba changad to 


Aj^(X^)y^^^+'i'’^(Xpua^)My|^4f’^(Xl) MyJ^^ - 0 . (4.21) 


Following tha aaaa llnaa aa tha baginning of tha proof and uaing tha ahlft- 
lng proparty tha following chain rula la obtalnad. 


dU(XJ 


Aj^(Xi)y^^^fV’^(Xi)U(Xj.)My^*"^^4-'i'^(Xi) My^^^V ... 


... + 't'^(Xj^) (1^2)1 ^ 

(4.22) 

Equation (4.14) of thla thaoran la obtalnad by racognlzlng that 


,(A+1) 


{ill) ^ w ~ ) 0*1 


thua tha thaorem la proved. 


Tha laft ganarallzad latent vectora aaaoclatad with aatlaflaa a 



alnllar chain rula that la given In tha next theoran. 
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Thtoria 4.2 L«t Aj^(X) ba daflnad aa in (2.1), & aat of laft latant vactora 

iotm a laft Jordan chain aaaoclatad with tha 
latant root and tha prlnary latant vactor . Tha chain rula if glvan 
by 

-1) X -<i-2) + 

^2‘ dX^ ^ 

*•*■’' (i-l)l dX<^X) 

(4.24) 

whare hj la tha langth of tha Jordan ch*.in. Tha vactora for Kfii<hj ara 

( 1 ) 

tha ganarallzad laft latant vactora of tha primary latant vactor z£ . 

Proof: Tha proof of thia thaoram followa along tha aama linaa aa thaoram 

4.1 whara tha chain of ganarallzad left alganvactora aa glvan in (3.6) la 
uaad with tha corraapondlng valuaa of Zq^^ . Tha alamanta of ara 
locatad in tha and antrlaa of 



Aa can ba aaan in tha proof of Thaoram 2.1, aoma intaraatlng ralatlona axlat 
batwaan tha gansralisad alganvactora of F and tha ganarallzad latant vactora 
of Aj^(X) . 


Connactlon 4.4 Lat y^^'' , . . . ,y^^ 

of F and lat y^^\y^^\. . . »y^ ^ ba tha 
latant vactora of Aj^(X) aaaoclatad with 
latant vactora ata ralatad by 


bo a aat of ganarallzad alganvactora 
corraapondlng aat of ganarallzad 
tha root X^:: Tha algan^<actora and 


. I, a 4. 

Foi - U(Xj^)Myj^ + 


fiV ^a-1) ^ 

dX “^1 ^ 




• * • * (A-1) I 

1 “ 1, 2,... ,hj^ , 


dX 


(4.SS) 
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Proof ; Tho oigonvictor !■ o fTjnctlon of woo provloualy atatod 

In eonnoctlon 4>2. Now ualng y^^^ aa glvan in (A. 19) and aubatltutlng y^^^ 




(4.26a) 


and ualng (4.20) 


v^2) . „a j. (1) 

yQi - U(Xi)^^r|_ + dX "yi • 


(4,26b) 


.(3) 


Tha ganaralizad elgenvactor y^^' can ba obtalnad from tha chain rula 
(F-X^l)y^^^ ■ y^^^ and following tha arguaanta aa In (4.19), with 




(4.27) 


( 2 ) 

SubatltUClng eha vector y^^^ aa glvan In (4.26b) in tha laat axpraaaion than 


- ua ^M,r^3) . ^”^^1^ „ (2) . (1) 

yoi - U(X^)My^ + —55;— My^ + 2f ^^2 %1 ‘ 


(4.28) 


Tha remaining alganvactora in tha chain are obtained from tha axtenalon of 
(4.28). WV 


Tha left ganaralizad alganvactora of F can be obtained from tha left 
generalized latent vactora of A^^(X) aa given In tha following connection. 

(2) (3) 0<j) 

Connection 4.5 Let generalized left 


01 


(hj 


(2) (3) i 

aiganvectora of F and let z£ be the corraapondlng aet 

of ganaralizad left latent vactora of Aj^(X) aaaoclatad with tha root \^. The 
left aiganvectora ted left latent vactora are related by 
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‘Oi 


^' V *1 * iX *1 ■ 


(4-1)1 


(4-1) 

3L.(1) 


^ * X|ii2f*tii 9^^ 


(A.29) 


Proof ; The left eigenvector as a function of waa pravloualy ob- 

T (2) (1 

talned in connection 4.3. The left eigenvector aatlaflea (F-X^I) 

(2) 

where the elemanta of the left latent vector z£ ^ are at the and (Vj^4V2> 
locatlona of z^^^ and thua by alaple aubatitutlon, z^^^ can be repreaented 
by 

*of " (^-30) 

where (ah U(X^))^ la the ahlftlng upwarda on the tranapoae of U(l). Thla 
ahlftlng matrix haa the property that 


1 

(4+1) 


(ah U(X)) 


T d 


(4) 


dX 


m 


'I'(X) 


(4.31) 


Ualng aa glvan in (4.12) and the above property with (4.30) » then z^^^ 
la glvan by 



'F(X^)z[^^ 


dnx,) 

+ ___ 2^ 


(4.32) 


Succaaalve application of the chain rule of the left generalized eigenvectora 
and (4.31) glvca equation (4.29) aa atated In thla ci.?nnectlon. WV 

It haa been ahown that the eigenvectora of the obaerver atata form can 
be defined In terma of the latent vactora of the row-reduced form Aj^(X) . 


V ■■ f ■ * 

Xt eh«r«for« followi that th« algravactora of tho oboorvoc atato fon can 
ba obtainad diractly froa tba knoirladta of tha lataat roots aad latant vac- 
tors of 


page is 
quality 


5, An Example 


omaiNAL 

OF POOR 


An lllustrmtlv* txinplft Is glvsn. Lst ths roir-rsducsd font 

of s Isabds Mtrlx that has boan subjactad to a unlaodular aMtrlx 


V^> 


X^-A^+7X+3 8X+5 


-3X-I-3 


X^-6X+5 


(5.1) 


and than ax|»rasslng tha abova rov-raducad laid>da aatrlx Into tha fon glvan 
In (2.1) 


V^> 



1 

>• 

0 

1 


1 0 


» • 
X^ X 1 0 0 

- 







1 

0 

>' 

10 

1 


0 1 
> ■ 


0 0 0 X 1_ 



0 

8 

5 

-6 

5 


(S.2) 


vharait can ba notad that ■ 3 and V 2 ■■ 2. Tha obaarvar canonical fon 

is 


F - 


110 0 0 
7 0 1 -8 0 

3 0 0 -5 0 . 

3 0 0 6 1 

3 0 0 -5 0 _ 


(5.3) 


Tha datanlnant of tha lanbda aatrlx la |Aj^(X') | * X(X-l)(X-2) , tharafora 
Xi 0, X 2 " 1 and X^ * 2. Tha prlaary right latant vactora Can ba cal- 
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eulat«d using (4,7s) with 


Aj^(0)y[^^ - 0 A^U)y2^^ - 0 " ° 

and tha latant vactors ara 



has two ganarallzad latant vactors that can ba calculated using (4.14) 


'3 ■ “ “■» A*«>T3 * ^3 

( 1 ) „ 

+ -r s — y» ■ 0 




(5.5) 


and tha ganarallzad latant vactors ara 



Tha right alganvactora of (5,. 3) can ba calculatad ualng (4.10) whara tha 
unlaodular aatrlx u(X) and M ara 


U(X) 


1 

0 

0 

0 

0 


1 

0 

X 

1 

0 

0 

0 


-1 

0 

CM 

r< 

X 

1 

0 

0 

M - 

7 

8 

0 

0 

0 

1 

0 


0 

1 

0 

0 

0 

X 

1 


-3 

-6 

_ 




> 


a 

. 


(5.6) 
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Th« •Igtnytcr.ori ar« than obtalnad wl^h 



^5 


"l3 ' 


• ■ 

1 

3 


0 


1 

^11 

„<1) 
^02 • 

11 

(1> . 
'03 

1 

3 


-10 


-1 

-3 

■ 


11 _ 


1 

b .a 


Th« |•ll•ralls•d alganvactorB for 



v«r« 


obUinad froa (4.25) 


- U(2)My^^> + 
and 

y<« - U(2),^<*> + ^ + 1 My<» 

dX 

( 2 ) 

whara y^^^^ and y^^' ara found to ba 


-1/3 


1 " 

2/3 


2/3 

0 

■'03 

1 

0 


-1 

0 

« 


1 . 


f 


Tha iaft latant vactora can ba calculatad ualng (4.7b) aa 


aJ(0)x^^^ - 0 . - 0 and Aj(2)a^^> - 0 (5.8) 


than 
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.( 1 ) 




1 

-I 


<1) ® 

*9 • 

L 1 


g(i) , 


[: 


and tha two gadarallsad laft latant yaejtora ara obtalnad from (4.24) 


A!(2)a^2) ^ (1) , 

V^^*3 ^ dX *3 




T M {« A ■ U » 


2 ^2 "3 


<5.») 


Tha laft ganarallzad latant vactof/i; ara 


.( 2 ) 


■ - 1/2 ■ 


- 3/4 ' 

_- ll /2 

3 

.- 7 / 12 . 


Tha laft alganvactora ara calculatad fro* (4.29) and thay ara 


01 


• « 
0 


m m 

0 


*4 

0 


0 


; 2 

1 

a<« - 

*02 

0 

- 

03 

1 

0 • 

1 


14 

.-1 . 


1 

m 


[7 j 


Tha laft ganarallzad alganvactora for a^jj^ ara than obtained from 


^(2) 

'03 


n2).f 


and 


- 


“03 




dX 


dX' 


(5.10) 
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with and 1^3^ livan by 



* 

2 


taaa 

0 

. 

*03 

«2 

- X /2 

- 3/4 

> 4 

-20/3 

-11/12 _ 


-7/12 


It !• obvious that tha aigsavaetora of tha obaarvar font can ba obtalnad 
froB tha latant vactors of tha roi^raduead form. 
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6. Concluiloni 

This ch«pt:«r has axplorad th« clganvalua-alganvactor ifciructur* of ob- 
sarvir canonical fora and Ifc has baan shown that this structure la ralatad 
to the simpler underlying structure of the corresponding lambda matrix. 

This work generalizes the prevloua results pertaining to a special class of 
lambda matrices In which the row Indices are taken to be equal. The special 
structure of the eigenvectors has potential for simplifying computations 
concerned with their evaluation. Dual rasulta also clearly exist for right 
fraction descriptions, column reduced lambda matrices and the controllmr 
canonlcel foraa. 


CHAPTER 4 


PARTIAL FRACTION EXPANSION METHODS FOR MATRICES 
OF STRICTLY PROPER RATIONAL FUNCTIONS 

Chapceri 2 and 3 daacribaa tha algabralc thaory of latent roota* vactora 
and projectoira, The latant projactora of the lambda matricaa are the matrix 
raalduaa of tha partial fraction expanalon of [A(X)]~^i aaa Section 1 of 
Chapter 2. Tha latent projactora are Important In obtaining the time domain 
iolutlona alnca tha Invaraa Laplace tranaform can be ased once the partial 
fraction expanalon la known. Tha latent projectors or toatrlx raalduaa can 
be computed by either using the latent vectors or by tha more claailcal 
procedure of computing tha raalduaa of the partial fraction expansion. 

Tha latter method la not reconnanded for a lambda matrix with hlgh-ordar 
matrix coefficients although algorithms for the procedure will be given 
here for completeness. 

This chapter extends several partial fraction expansion methods for 
matrices of strictly proper rational functions. These methods are 
1) eigenproj actor, 2) Chen and Leung and 3) escalation which ate first given 
for the scalar problem and then extended to the matrix case by using the 
Kronacker product. The three methods are competitive for high order functions 
although the escalation method Is the most efficient for scalar functions. 
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1. Imtroductiott 

It Is frtqusntly nscssssry in systsa snnlyiis to find ths psrtiil 
fraction expansion of a aatrix of proper rational functions [1] and [13j. 
Several Mthods for the expansion of scalar rational functions have been 
dafcribed in racantly published papers* The paper by Kami and Etter[17] 
was aaong the first to show that derivatives of the functions ware not 

I 

necessary in the evaluation of the residues of the expansion* Chen and 
Leung [IS], es well as Mahoney [19]» also avoid derivatives and describe 
algorithms that are quite general. The paper by Shahzadi[20] described a 
procedure that eliminates the necessity for the differentiation of rational 
functions but derivetivaa of the numerator and denomlnatOs' must still be 
taken separately. The authors have developed an algorithsi which is quite 
similar to the procedure given by Chen and Leung, The algorithm developed 
by the authors will be described In this paper and will be denoted as the 
eigenproj actor method [ 113 . 

The procedure given by Mahoney la the most efficient in a numerical 
sense for scalar functions but the method does not clearly reveal the alge- 
braic structure of the partial fraction expansion procedure. The Chen and 
Leung method, as well as the eigenproj actor procedure, is basev: an the 
algebraic structure of the partial fraction expansion. The Kami and Etter 
method is algebraic but the referenced version had a limitation in that the 
order of the numerator polynomial should not excaed the order of the product 
of the distinct roots of the denominator of the rational function. 

This paper will first give a brief mathematical description of the 
three asthods 1) eigenproj ector, 2) Chen and Leung and 3) escalation 
(Mahoney). The elgenprojector msthod has not been described in the liter*- 


67 


original page is 

OF POOR QUALITY 


atur* and for that raaaon, tha davalopaaint will ba nora axtanaiva than for 
tha othar two procaduraa. Tha Chan and Laung nathod will ba davalopad and 
tha Maltonay procadura will ba auanarlzed. In Sactlon 2^ tha daacrlptlon 
of tha Mthoda will ba liodtad to acalar functiona. 

In Sactlon 3, tha nathoda ara aztandad to Batrlcaa of propar rational 
functiona of the fora 


H(a) 


C,a“"^+C,c“”^+... 

1 z n 

n, j n—1, , . 

a +d,a +. . . • • t-rd 

1 n 



( 1 . 1 ) 


whara and tha d^^ coafflclanta ara acalara. Tha partial fraction 

axpanalon will have tha fora 


H(a) - I I 
1-1 1-0 


*‘ 1.1 

(a-X^) 


m: 


( 1 . 2 ) 


whore la tha multiplicity of the root and r la the number 

of dlatlnct roota of dCa) . Zadeh and Deaoer II] diacuaa tha axpanalon of 

functiona of the type given above, but the determination of the raaiduea 

K. 0 In their work followa the derivative procedure. 

*•* 

An example will be given In Section 3 to llluatrate the three methoda 
aa well aa the number of operatlona required where the usual aasumptlons 

3 

that the Inveralon of an nxn matrix requlrea n operatlona, the aolutlon 

3 

of n equatlona requires n /3 operations [21] and multiplication of two nxn 

3 

matrices requires n operations where an operation is defined to ba one 
scalar multiplication together with one scalar addition. 
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2. KathQn&tlcal Oavclopnant fur Scalar Casa 

■) 

a) Eiganprojactor Hathcd 

A atrlctly propar rational scalar function H(a) can always ba diarac- 
tarlzad by tha trlplat (A,B,C). If H(s) la glvan by 


1 

H(a) " 


C,S •►CjS +. -fc^ 

n. . n-1. ^ " dCa)” 

i. n 


( 2 . 1 ) 


than H(s) can also ba daflnad by 


H(a) - ClaI-A]’^B 


( 2 . 2 ) 


whara AeR°^^» and CeR^’^". Tha atructuraa of A, B and C ara glvan by 

I9J, 


A ■ 


b". 


and 


-d 


n-1 


0 9 ... 


1 . . ‘ • 


0 

-d. 


• 0 
• 1 
-d. 


(2.3) 


0 , 1 ] 


(2.4) 


®n-i» ••• •••» ®2* ®1^ 


(2.5) 


whara the poles of H(a) are equal to the eigenvalues of A. 

Zadeh and Desoer have shown that H(s) as given In (2.2) can be repre~ 
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■•ntad by tha axpanalon 


r “r^ P 


H(a) - C{ j; I 


l.A 


i»l i-o (a-X^) 


£+1 


>B - I I 


r k 


lA 


1-1 1-0 (i-xp 


1+1 


( 2 . 6 ) 


vhara tha nxn matrlcaa ^ ara tha Matrix raalduaa or alganpro j actors » and 
l^ha scalar raslduas. It Is assumad that H(s) has r distinct polas 


and aach pole Is rapaatad tlnas such that 


ft •; 

M 


i* «4 - ti . 


(2.7) 


1-1 


I ' 


Tha alganproj actors, can ba detaralned from [sI-Aj"^ idilch can ba 

rapresentad by tha following similarity transformation 


taI-Al“^ - M[sI-J]"V^ 


( 2 . 8 ) 




I f 
i I 


where M Is tha eigenvector matrix of A, and J la tha Jordan canonical form. 
Tha alganproj actors can ba defined In tha usual sense as 

1 • 0,1, ,.^-l (2.9) 

-1 

where the structure Of [sI-J] la , [10] ^ 


,-l 


[sI-J]'-" - Block Dlag{[al-Jj^]"^,[sl-J 2 ]“'*' .[sI-J^]"-^} (2.10) 


-1 


,-li 


It is not difficult to show that the Inverse of sl-J^ Is given by 
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-1 




♦ 

0 


<-^l>“'. 


(i-X^) 


-I 


(■-X^) 


-1 


•’ 0 




(.-X,) 

(•-x^) 


-2 


. (2.U) 


-I 


It follows from (2.9) that tha alganproj actor P. , will hava tha fol- 
loving form: 



. . I ri 


- M(ah F. . ,Im"^ (2.12) 
x»m^“x 


whara the m^^m^ Interior block has the same location as the 1th Jordan 

block In J. The matrices [sh F. .] will be defined as the shifting 

i,m^-x 

matrices. The next elgenprojector, F. „ for the repeated eigenvalue X. 

x,m^-4 1 

Is obtained by shifting the ones to the next lower diagonal of the Jordan 
block. The remaining elgenprojectora for the repeated eigenvalue X^ are 
generated by shifting the ones along the remaining lower diagonals until 
the main diagonal Is reached. If the eigenvalue Is not repeated, then a 
one Is placed on the diagonal position corresponding to the Jordan block 
for chat eigenvalue. 
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0 

Th« «lg«nv«ctor matrix M of A !■ conatruetod from the roota of d(a) aa 


“ 




1 i 
1 ) 


1 

1 II ^ 

0- . ‘ 

... 0 

1 


1 1 i»i 
1 

1 . 

• 

: 

1 ’ 
1 1 

s 


ills 

2X, 

• . 

’ 1 

II 
1 1 


• 

• 

(1 : 

• 

• 

* 1 

1 1 
1 1 

9 

j^n-1 

[h 

" x“^ 

MX, 

("-Sx;-* 


1 1 
i 1 



(2.13) 


where the eigenvector for a non-repeated root la aa given on the first 
column of <2. 13). Since is a pole of H(a), the eigenvector matrix M can 
bo constructed directly from the poles of H(s) . 

The olgenprojectors given in (2.12) are related to the scalat residues 
of the partial fraction expansion in (2.6) and the scalar residues are 
obtained by premultiplying the elgenproj actors by vector C and postmulti- 
plylng by the vector B as given in C2.4) and (2.5). The raaiduas are than 
given by 


'‘l.l • “l.l‘ ■ “f'" 

and the only difference between residues for different eigenvalues is in 
the shifting matrix [sh which is a matrix of onea and zeros which 

selects the proper columns of CM and the proper rows of m'^B. 

The computational procedure for the elgenproj actor method ia detarmln- 
Ining the residues of the partial fraction expansion of HCs) is as follows: 
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(I) S«t up th« •Iguayuctor Mtrlx M, froB knowludg* of th« roots of 
d(g) M givon in (2.13). 

(II) Lot CM B ond M~^B ■ and coBputo x and y from C and B 
no glvsn In (2.4) and (2.5). 


(Ill) Sot flags regarding the aultlpllclty of each root to carry out 
the shifting proceos on [ah 

(Iv) Compute ^ from ^ - x[sh ^Jy. 

n3 2 

The algorithm will require approximately y- + 2n - n operations for 

distinct roots provided that the flags set In (Iv) are used. For a single 

n3 2 

root with multiplicity n, the count goes to approximately •¥ 3n . If n 

3 

Is reasonably large, the n /3 term will dominate, thus this value Is on 
approximate count. 

The following example is given to show the procedure of the elgorlthB. 
Let H(s) be given by 


H(s) 


3s>854-9 _ "‘10 . . 


20 ^ 21 


a^+7s^-fl5s+9 e+1 s+3 (a+3)^ 


(2.15) 


for which the eigenvector matrix M con be reconstructed using (2.13) with 



1 1 o" 


11 0 

M - 


- 

-1 -3 1 


x] 2 X 2 ^ 


1 9 -6 


The X row vector and the y column vector ore now computed where 
C - [9 8 3] end B*^ - [0 0 1] with 
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X - CM •• [9 8 3] 


1 1 
-I -3 
I 9 


0 

1 



[4 12 - 10 ] 


•nd 


1 1 0 

-1 

0 


"1/4" 

-1 -3 1 


0 

m 

-1/4 

1 9 -6 


1 

vm wm 


-1/2 


and for thla axaqplo tha ah utrlcaa ara 



1 0 

0 


000 


0 

0 0 

■“ '1,0 ■ 

0 0 
0 0 

0 

0 

*•' '2.0 ■ 

010 

001 

'2,1 • 

0 

0 

0 1 
0 0 


In programalng the algorithm for a digital computer » the flaga would be aet 
to Indicate which combination of elementa In x and y are neceaaary to cal- 
culate the realdue. The partial fraction expanalon of (2.15) la 


H(a) 


1.2 6 

+ 2 

8+1 8+3 (a+3)‘ 


which la the desired expanalon. 


(2.16) 


b) Chen and Letmg Method 

The Chen and Leung algorithm, although similar to the elgenprojector 
method, differs In that all the realduea of the partial fraction expansion 

are computed slmultaneoualy. 

1 : 
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Utilizing (2<1) znd (2.6), It follows that 


1 2 n (s-X^) 


r “r^ K.2 


(s"+d^t“’V . . ,+d^) (2.17) 


Tha ratio of tha two polynonlala on tha right aids of (2.17) can ba wrlttan 
In tha ganaral fora 

Substituting Into (2.17) snd equating tsnu In s^ glvas 

'2 ■ ■*■ Jj ’'l.l * 

or In general 

" Ji'^l.O ill'll.! + ..• 


( 2.10 


(2.19) 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


where 


d^_j^(X) - X^"^+dj^ X*^"^+ 


If a pole is not repeated, then k. - 0 for all 1? i. 

XfjL — 

Assuming that all roots of d(s) ■ 0 are distinct, then (2.19)->(2.22) 


gives the matrix eqtiatlon. 
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‘‘1,0 


1 1 1 1 ) 1 “ 

X X ^ • A 

^1 -2 III \ 

1*2.0 

a 


« 

• 

■1 

: 1 1 r 

a a 

« 


• 1 1 1 • 

‘‘n,0 


xrM 1 1 

X A n 
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■“l 

0.. . 

a a a a 

0 

-1 



a 


a 



"i 

i 

a 

• 

a 

a 

a 


*^2 

• ♦ 

a 

a 

• 


a 

• 

a 

a 

a 

4 

a 

' . 
a 

. • 

1 

a 

f 

0 


♦ 

a 

a 

Jn-l 

a 0 . e 

. . ’ dj 

X 


c 


(2.23) 

The netrlx equeclon of (2.23) la the daalrcd exprasaion for conputlng the 
raalduoa. If the roota of d(a) art rapaatad, tha alganvactor natrlx M, 

aa ahown on tha right alda of (2.23) will taka tha ganaral fom given In 
(2.13) with tha raaldua vector on tha left side having tha raalduaa ^ 


In aaquentlal order with , 1 « . . . , 1-1 . 

Tha cooputatlonal procedure for the Chan and Leung algorltha la quite 
ainpla and atraightforward, the algorithm la: 


I) Set up the eigenvector aatrlx M and the coefficient aatrix D. 

II) Find the aolutlon of the aet Of equatlona froe D and the coaffl- 
clenta vector. 

III) Find the aolutlon for the aet of equatlona between M and the 
vector found In (11) . 


The prevlouB exaaple will be conalderad ualng thla method. The eigen- 
vector matrix M la the aame aa before therefore 


‘‘lO 


”l 

1 

0 

-1 

1 

0 

(T 

-1 

‘‘l 

‘‘20 

- 


^2 

1 


‘‘l 

1 

0 


*=2 



A 

^2 

2X^ 


_-‘2 


1 


J3_ 
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and upon subsicutlng valuaa Into th« «quaeion» tho raalduaa ara 


*‘io 


1 

1 

0 

-1 

1 

0 

0 

-1 

’ 3 ” 


1 

*^20 

- 

-1 

-3 

1 


7 

1 

0 


8 

- 

2 



1 

9 

-6 


15 

7 

1 


9 


-6 


idklch are correct. 

The number o£ operatlona required for computing the rcalduea ualng 

2n^ 2 

this method ia approximately n -2n for all diatlnct roota and 

3 

2n*^ 3 2 

7 n approximately for a aingle root with multiplicity n, ao for 

^ 3 

2n'^ 

large valuea of n the term '-j- will dominate, thua thia value ia an ap- 
proximate count. 

c) The Eacalation Method 

Mahoney haa recently publiahed a paper on partial fraction expanaions 
by the escalation method. The method is quite efficient although the alge- 
braic structure of the process is not as apparent as in the previous two 
methods. The mathematical theory will not be given in this paper and only 
the algorithm will be presented with the extension to the rcalduea of a 
matrix of proper rational functions in the next section. 

Consider a scalar function as given in (2.1) where c^ and d^ are 
scalars. The basis of the method proceeds as follows; use Homer's algor- 
ithm [23] to expand the numerator of H(s) into the form 

n-1 i 

N(s) - 6. + j; 6. n (s-X.) C2.24) 

” i-1 j-1 J 


and, therefore, H(s) can be expressed as 
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0 + 0. n (.-X ) 

«(.) - ».(.> - — — - 

J-1 ^ 

wh«r« th« abovs tiqiraasion can ba dacoqpoaad Into 


" (a-Xjj) 1:=^ 


A - 1. 


<2.26) 


Tha function vill have a aladlar atructura aa (2.25). Lattlng v ba 

an indax on tha roota of d(a)« than tha partial fraction axpanalon (2.26) 


can ba vrlttan aa 


^(4) ,(« , (4) 

H,(.) . r + 1 -“-T 

^ j -1 (8-^x.)^ i-1 J-1 i»-Kr 

^ l|‘v ^ 


a. aaj 

+ I I" 


l<A<n 


(2.27) 


tdiara and la tha nultlpllclty of tha root X^ In H^(a). Aa 

can ba aaan In (2.26) whan ia la allalnatad the multiplicity of 

that root In Hj^^^(a) la reduced by one, ao tha acalatlon method la baaed on 
that principle. If X^ la taken out In (2.26) than tha aultlpllcltlaa of 
tha roota will ba aa 




a ■ X. 


KV<r 


(2.28) 


, (A-1) . , (A) 

h \ 


1 - 1,2, ,m. 


1 f* V . 


Tha proceaa tarmlnatea whan A-n and - m^, 1 - l,2,...,r. The coaffl- 


(AV fA^ 

clanta and are computad from 


,(A-1) 

b(a) . !ii_ 

A -X 


j - L^^^>1 


C2.29) 


. • ssi. 
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and 


M) 


'ij 


-(JUl) .<l) 


h”K 


j • • • • »2|1 


(2.301 


which Mahonay ealla tha allan coafflclanta. Tha coafflclanta glvan by 


,<A) 


ii) 


'VI 


^1-X “ ®1X 

iFr 


(2r^l 


.(A) . .(A-X) 


'vj 


’V,J-X 


J ■ 2y3|*..»L 


(A) 


(2.32) 


ara danotad aa cha naClva coafflclanta. 

The ascaXaClon procadura la afflclant but tha pcograaaiing affect to 
IppXanant tha aXgorlthai la graatar than althac of tha two pravloua aXgoc- 
Ithaui. Tha conputatlonaX procadura la aa foXXowa; 


1) Uaa tha nuaMcator coafflclanta and tha Hornar'a algorithm to 
coaqmta the coafflclanta 0^ In (2.24). 


11) A tableau la then fotnad from Equatlona (2.29)~(2.32) and tha 
raalduaa will be obttiilnad In tha Xaat row of tha nxn tableau., 




The exaiaple given before will now be conaldered ualng thla method. 


Tha Hornar’a algorithm la uaad to axpraaa tha numerator In tha form ahown 


In (2.24). For n(a) - 3a^*f8a-l'9, the Hornar'a achama glvaa 


3 

3 

3 

3 


3 

5 

-4 


■ -1 


Xj - -3 


■ i 

■ '> i 

; I 


r 

I . 
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1 

thus I n(f) •• 4-A(e4’l)4*3(a<fX)<B'^3). Th« tabltau la nov conaCructad ualn| 
thaaa valuaa and Eqa. (2. 29)- <2. 32) with tha tableau glvan aa 



*11 

®21 

*22 

r-1 

A 



r-2 

2 

-6 


r-3 

1 

2 

•' 1 ' 


*^10 

o 

CM 

*‘21 



whara tha raalduaa of tha axpansion ara found In tha laat row with agroa- 
aant to tha valuaa glvan aarllar* 

Mahonay haa aatlnatad tha nuabar of oparatlona for tha procadura aa 
approximately Thla mathod la therafora more efflclant whan axacu- 

tlpn time la tha only factor conaldarad for digital Inplamantatlcn. Tha 
procadura la claarly auparlor to tha other mathoda whan manual calculatlona 


ara mada. 
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.1. Extanslon of tht Methods to. Matrices of Proper Hationel Functions 

The three procedures uhlch have been cXasslfled as^ 1) the eigenpro<> 
jector« 2) the Chen end Leung end 3) the esceletlon aethod will be extended 
to strictly proper aetrlx function defined es 


jj; H(s) - 0 . 


(3.1) 


It will be essuasd that the function has the £orm 


H($) 


* BL 

8“-fd,s“"^+ +d„ 

1 n 


(3.2) 


where with scalar coefficients in the dsnonlnator. The Kronecker 

product will be used to extend the aethods as well as to sdnialse the 
arithaatlc operations required for each SMthod. 


a) Eigenproj actor Method 

It can be shown [24] that a rational matrix functions as given in (3.2) 
can be expressed in the fora 


H(s) - C(Is-A)"^B 


A - 


•d I -d -d^ 

n q n-1 q n-2 q 


(3.3) 




• 0 


-“iS 


- A0I. 


(3.4) 
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■ I0q*0q»*»* ’••*®q*^q^ " 

and 

(3.5) 

C ■ •••*^2*^1^ * 

(3,6) 

Using the definition of Kronecker products (25] and [26], 
in (3.4) can be defined as 

the nqxnq matrix 

I - (M© Iq) ( J © Iq) © Iq) 

where J and M are as previously defined 

(3.7) 

It follows directly from (3. 3)- (3. 7) that the matrix of rational func^ 
tlons H(s) can then be expressed in the form 

H(s) - C(M© Iq)((sl-J)”^© Iq) vM"^0Iq) (B© Iq) 
or 

(3.8) 

H(s) - C(M© Iq) ((sI-J)’^© Iq) (M"^B0 Iq) . 

(3.9) 

The matrix residues of the expansion in (3.2) can be 
from (3.9) as follows 

easily obtained 

*^i.i “ G(W©Iq)'(*h j^®Iq)*(M‘H®Iq) 

(3.10) 


wharc [sh F. «] la as defined In (2.12). The cos;>utatlonal procedure Is 

l,x> 

then as given in the scalar except x ■ C(M©I^) and y ■ M where x 
Is a pxnq natrix and ^ is an nqxq matrix. Note also that in this case the 
last term in (3.10) can be cooiblned with the shift operation and thus 
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‘‘l.i - =(«©!,) <•'> 

To llluBtratOi conoidor tht oxa^plo 


3 7 

2 


8 

24 

15 

a+ 

9 

17 

9 

—6 5 

10 


-31_ 

82 

45, 


-45 

51 

47 


•‘+7i‘-fl5«+9 


whoro th« root! of tho dononinotor ar« ■ -1 with * X X 2 
1 B 2 ■ 2. Tho eigenvector matrix M ia given by 



*"1 1 o"" 


1 1 o' 

M - 

Ai Xj 1 

■ 

-1 -3 1 


fl «2. 


19-6 


and CCmQi^) ie 


G(M©l3) 


4 0-4 

-20 24 12 


I 12 8 -18 
I -6 0 2 



B will be given by -• [0 0 1], therefore 


(m“V©I3 ■ 


Now using the expreaeion for the reaiduea it follows that 


1/4 

-1/4 

■ 1/2 


© 1 , 


(3.11) 

(3.12) 
- -3 with 
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^10 “ " 

K 20 ■ *l»h '20®^3^^ “ 

iCj_^ - 5(.h FjjQIj]? F‘ ■* 



Again* there la no need to construct the [sh aatrlz as the daslrad 

operations ara carried out by using tha flags. Tha matrix partial fraction 
expansion of (3.12) la then 


H(s) 



(3.13) 


which Is correct. 

The number of arithmetic multiplies and divisions Is approximately 

3 2 9 9 3 2 

+ 2pqn^4n* for a low count and + 2pqn + n for tha high count. 

b) Chen and Leung Method 

The derivation for the scalar case carries over to the matrix of proper 
rational functions In on analogous manner to the elganprojector case. Eq- 
uation (2.23) Is modified by performing a Kronackar product of M and D with 

I • thus for the distinct caaa 
P • 


'• - 4. 
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*^1*0 


*^2*0 

• 

f 

• 

- 

* 

* 

n *0 



^2lp 


^1 V ^2 



-1 

I 


1 ** 


{111 


II n p 

1 1 • 

• 


1 1 • 
i 1 ' 




M n p 



d.I I. 

IP P 


d , X • • • 
n-1 p 


“p 

d,I I 

1 p p 


I I 1 


(3.14) 


vh«r« aad I and 0 ar« pxp. Using ths properties for Kronecker 

products ss given in the Appendix* theri (3.14) csn be rewritten ss 




V 

■'2.0 

» 

- ( M ‘ V ^© Ip ) 

<=2 

« 

e 

1 

n *0 


s 

c 

L“J 


(3.15) 


where M end D ere identlcsl to the mstrlces for the scslsr esse. 

Using the sssm exsople ss before* the mstrlces M snd D sre for the 
example 



1 1 0 

-1 

o 

r-i 

L _ 

-1 

1 -1 1 


-1 -3 1 


7 10 

1 

" 4 

3 1 -1 


19-6 


15 7 0 

mamm «k>J 


-18 6 


snd using (3.15) the following expansion Is obtained 


*^ 1.0 


1 0-1010 

0 10-101 


3 7 2 

- 6 5 10 

*^ 2*0 

1 

4 

3 0 10-10 

0 3 0 1 0 -1 


8 24 15 
-31 32 45 

" 2*1 
L ^ 


-18 0 6 0 -2 0 

0 -18 0 6 0 -2 


9 17 9 

_-45 51 4^ 
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whttrtt th« rialdu«f found from th« oparatlona ara tha aaaa aa glvan (3. 13). 

3 2 2 

The method requirea about 2n 4pqn for a low count and about 
3 2 3 2 

2n 4pqn ^ ^ which are ulakoit tha aama for largo matricaa. 
c) Eacalatlon Method 

The laat method to be conaldered for axtenalon to the expanalon of a 
matrix of proper rational functlona la the eacalatlon method. The uaual 
procedure would be to extend tha acalar caaa but where matricaa era uaad 
In tha procedure aa followa. The Homer’ a algorithm la uaed to caleulete 
the 0^ uatrlcea of 

n-1 1 

N(a) - $0 + I ^ <3.16) 

0 iil ^ J-1 3 


where 0^, . . . Equ&tlona (2.29)-(2.32) would thfu be uaed In a 
matrix form of the tableau. 

For the example preaented In (3.12), the Homer 'a algorithm glvea 


N(a) 


T4 0 - 4 “] ^ r.4 -4 fl f 3 7 2 “] 

-20 24 12 \ -7 12 5 -6 5 1^ 


<a+l)(a+3). 


Conatruct the tableau In the aama manner aa the acalar case 


given 


11 


,(r) 

*21 


|(r) 

*22 
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th« iMt matrix raw of tha tablaau glvaa tha raslduaa and thay agree with 
C3.13). 

The arithmetic count, conaidering only multiplies and divides, is 
2 

approximately pq(n 4ti/2). The method is more efficient than the other two 

!i 

but the efficiency tends to be less drastic than that for the scalar case 
particularly for p and q approximately equal to n for large n. 

The Mahoney procedure can be modified so that fewer operations are 
required when compared to the matrix formulation given above. Consider 
the matrix of rational functions given in (3.2) which can be expressed In 
the form 

[c. ,0- ...... ] 

H(.) - „ ^ -s— 

•**+d,.** ••y.+d 

X n 

The functions 

sJ 

h.(s) - — J - 0,1 ,n-l (3.18) 

J s”>d,8“ +d 

X Q 


.n-l 

n-2 


©I, 


(3.17) 


appear in (3.17) thun knowledge of the n-expsnsions of (3.18) will provide 
the basis for the expansion of (3.17). The Mahoney procedure can be used 
to expand h^.j^Cs) and the tableau can be continued to find the expansions 
of * 

If h ,(s) can be represented by 
n**x 


»‘n-l<»> 


r 

I I 

i-1 £-0 


*^i.£ 


£+1 


(3.19) 


it follows that the expansion 
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“ j i-i 


T 

£-0 


,(j) 




A+1 


J ■ 1.2,, 


will contain coafflclants glvan by 


1,(J> 

^(j+1) , 




and 


(J+1) 

1 . 


A A 


^1 


ilJ±l 


A - m ^-2 .1 
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(3.20) 


(3.21) 


(3.22) 


which are the same equations given in (2.29) and (2^30) except X^*0. The 
latter two equations are then used to continue the tableau such that the n 
expansions of (3.18) are given In the last n rows of the tableau. 

Now It can be seen from (3.17) that the matrix residues mr. 

then obtained 


*^1.1 * ^^^1*^2' 


,Cnl 


( 1 ) 

!,«- 

( 2 ) 


v<“> 


©I, 


(3.23) 


where the vectors [k^ arc given In the corresponding columns of the last 
n-rows of the tableau. 

Considering the example In (3.12), the functions given by 


hj(-) - 


sJ 

(s+1) (e+3)' 


J - 0,1,2 
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■uat b« «xp«i<l«d. Starting with J-2, tha Bomar'a algoritha glvaa tha tahla 



and tha usual Mahcnay tableau for h 2 (a) la glvan by tha flrat thraa rowa of 
the tableau with the last two rowa giving h^(a) and hQ(s). 


r-1 

r-2 

r»3 


»(r) 

hi 

«<r) 

hi 

«(r) 

hz 

1 


-1/2 

-9/2 


1/4 

3/4 

-9/2 

-1/4 

1/4 

-3/2 

1/4 

-1/4 

-1/2 

j,<d) 

*10 

k(d) 

*20 

u(J) 

*21 


6, -.4 

6j-l 


The corresponding vectors of the expansion of C3.18) are 


p(l) 


1 


*2,0 


3 


7(1)“ 

*2,1 


-9 

*1,0 

4 

-1 

8 

*2,0 


1 

and 

k<2) 

’‘2.1 

1 

" 2 

-3 

k(3) 

_i,o_ 


1 


*2,0 
mm M 


-1 


k(3) 
1^2. ij 


••1 


and using these vectors In C3v23) will give the correct matrix residues. 


The operation coisit for the above modified procedure Is approximately 

2 2 

pqn -f2n -n/2 which Is less than the count for tha Mahoney matrix formulation 
as well as the other procedures, 
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4. Conclusion 

Three methods for computing residues of proper rational functions have 
been presented and they all have the property that derivatives for Che case 
of repeated roots are not necessary. The methods have been compared in 
operation counts and the escalation method is the most efficient for the 
scalar case. 

The eigenprojector and the Chen and Leung methods are based on the 
algebraic structure of the expansions which is elegant and easy to follow 
and the method can easily be extended to the matrix case by continuing 
the tableau and then by applying the Kronecker products. This method is 
shown to be the most efficient in operation counts. VRien the algorithms 
are implemented in a digital computer, programming attention should be 
given to the Kronecker products to minimize the number of arithmetic 
operations. The Kronecker products have been used in the development as 
a mathematical Cool but that does not mean that the equations are to be 
Implemented as given in the paper. 


CHAPTER 5 


ALGEBRAIC THEORY OF A( A ) - M + C A + K 

Chapters 2 and 3 of this report gave the mathematical development of 

the algebraic theory of lambda matrices. Lambda matrices of general order 

were considered and the properties of the latent roots, vectors and projectors 

were considered. The finite-element model of a structure will generally be 

2 

defined by a second-order lambda matrix Ml + CA-f K where M is the mass 
matrix, C is the damping matrix and K is the stiffness matrix. The latent 
roots of such a lambda matrix will generally be distinct except for the 
free-free modes which will be repeated and located at the origin of the 
elgenvaltie space of the associated state-variable matrix. 

This chapter discusses the properties associated with the specific 
lambda matrix that arises in the finite-element model, l.e., the second-order 
lambda matrix. The displacement, damping and Identification of the c ‘'.ructure 
is of interest in the drvelopment. The objective in this chapter is to 
consider the second-order lambda matrix utilizing the algebraic theory 
developed in the previous chapters. The spectral factorization of the second- 
order lambda matrix will be described as well as properties associated with 


the spectral matrices. 
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1. Introduction 

The previous chapter of this report gave the general algebraic theory 
of latent roots, vectors and projectors of laabda aatrices« The work did 
not consider the special esse of a second-order laabda matrix such as the 
one occurring in the flnlte-eleasnt model for structural analysis. The 
dynanics of a vibrating structure can be characterised by the second-order ' 
differential equation 


M 


d^x 


dt 



+ Kx - P(t) 


( 1 . 1 ) 


where is the mass matrix, CeR®^® is the damping matrix, IteR®*^® is 

the stiffness matrix, x(t)eR®^^ is the displacement ^ctor end F(t)eR®^^ 
is the force acting on the structure. It will be assumed that M, 1c, and K 
are symmetric positive definite matrices. 

• I 

The homogeneous equation 
2 

M ^+ C ^+&-0 ( 1 . 2 ) 

t 

can be transfonaitd to a lambda matrix by assuming that x(t) " Xq exp(Xt) 
or by taking the Laplace transform of (1.2) with zero initial conditions. 

The laabda matrix of Interest in the following work will be of the fora . 

A(X) - + CX + K . (1.3) 

This lambda matrix will have 2m latent roots that occur in complex conjugate 
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pairs, 2tt right latant vactora In conjugata pairs as wall as 2a latant pro- 
J actors that ara also In conjugata pairs. Tha noroallsad laabda astrlx 


xa) - + ox + K 


<1.4) 


will be takan as tha canonic fom a Inca tha assoclatad block coapoinlon 
■atrlx 


\ - 


0 





(1.5) 


la consldarad to be In canonic forn. As given previously, the algenvaluas 

of A are equal to the latent roots of A(l) and the elganvactors of A con- 
c c» 

tain tha latant vectors of A(X) as sub vectors. Tha latent projectors of 

A(X) are siiboatrlces of the elgenprojectors of A . 

c 

The algebraic development In tula chapter will be directed toward 
establishing the theory and structure of the latent roots. Vectors and pro- 
jectors of A(X) as given In (1.3) and the structure of C for either speci- 
fied mode dampings, or the location of the latent roots of A(X). The first 
step In the development will be that of transforming A(X) Into the canonic 
forn and using the canonic form for the algebraic theory of lambda matrices. 
Tha damped and undamped lambda matrices will than be studied and the alge- 
braic structure of each fom determined. 
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2. G«n«ral Xhaory of Laabda Matrlcao for Unduq;>«d Vibrating Syataaa 
The ganaral theory of vibrating ayataaa nnd laabda utricaa haa bean 
covardd by tancaatari [3]. Soaa of that nuitarial ia proaantad hare to pro- 
vide the proper setting for the devalopaanc that followaii 

The undanpad systaa that has no viscous or hyataratic daaping and the 
free vibrating syataa can ba characterised by the ■ aacond-ordar differen- 
tial aquations 


M + Kx ■ 0 
dt^ 


( 2 . 1 ) 


where x(c)er"^^. The syataa is undaaped and ths displacesanta x^Ct) will 
be sinusoidal therefore if it ia assuned that 


x(t) 



( 2 . 2 ) 


and equation (2.1) takas the fora 

(rMu^+K)q - 0 (2.3) 

2 

If y - - 0 ) , (2.3) can be rewritten aa 

(MM-HC)q ■ 0 (2.4) 

* 

where is a real, regular and synnetric aatrix pencil, [3]. A pencil 
of natricea ia regular if a) M and K are square aatrices and b) M ia nou- 
aingular. 
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K slapl* Mtrlx p«nciX of ordov m hao th« propactiao [3|. 

a) It is a ^regular pancil 

b) It haa n llnaarly Indapandant right and laft latant vaetora, 
and raapactlvaly. 

Savaral thaoraaa and corollarlaa froa Lancaatar [3] will now ba givan wlth^ 
out proofa. 

f 

Thaoraa 2.1 Latant vaetora corraapondlng to distinct latant roots of a 
ragular natriJE pancil ara linaarly indapandant. 

Corollary 2# 2 A matrix pencil of order u having a distinct latant roots is 
a slaple pancil. 

Thaoraa 2.3 A ragular aatrix pancil Mp-fK is aiapla if and only if for 
avary latant root tha aatrix HP'fK haa daganaracy aqual to tha aultlpllc- 
ity of p^. 

Thaoraa 2.4 If H and K ara raal syaaatric aatricaa and H is positiva da» 
finite, then Np+K is a siople matrix pancil. 

Corollary 2.5 Under tha assumptions of Thaoraa (.3.4) all latant roots and 
latant vectors of tha pencil Mp4-K are raal. 

Thaorem 2.5 The regular pencil Mp+K Is dafactlva if and only if there 

T 

exists a latent root p^ with a right latent vector y^ such that 
for all latent vectors of p^. 

Corollary 2.6 As a result of the theorems, tha matrix pencil Mp-t-K is thare> 
fora a ilmpla matrix pancil and has m latant roots with m linaarly Inde- 
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:;^ti n(lent latent vecfcors. 
The lambda matrix 


X(X) - + K . 


(2.5) 


will therefore have 2m latent roots X^ - The latent vectors 

for (2.5) will be real since (MX^+K)£R“^® and the latent vector for +Jw^o 
will be equal to the latent vector of All latent roots appear in 

imaginary pairs, there will be only m linearly Independent right 

and left latent vectors. Furthermore, since A(X) is symmetric, the right 
and the left latent vectors must be equal for a latent root X^, l.e. 



(2.6) 


where A(X^)’ j|^ ■ 0 and aF(X^)z^ ■ 0. 

The latent projectors for on undamped vibrating sy/^tem for a latent 
root Xj^ are given by 


- -T 


- -T 

h ^i 


i,0 ^ dA(X^) _ ^ dA(Xj^^) ^ 

*i^dX 


^i ^1 


dX ^1 


but since dA(X^)/dX ■ 2M X^, (2.7) becomes 


N - 


(2.7) 


( 2 . 8 ) 


Theorem 2.7 The latent projectors for an imaginary pair of latent roots 
of the lambda matrix for an undamped vibrating system occur in imaginary 
pairs 


\\ 
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(2.9) 


whace it has been assumed that the latent roots are arranged In the order 

♦ 

^21-1 " ^^10 ^21 " ”*^^10’ latent projectors for each pair are 

linearly dependent. 

The proof of Theorem 2,6 follows directly from (2.8). If A 2 i^l ” ^^10 
^ 2 i " “^‘**10 


^21 ,0 




1 /I ’U 


^^“10 


p * 

^ 21 - 1,0 


“^ 21 - 1,0 


( 2 . 10 ) 


Theorem 2.8 The latent projectors for different pairs of latent roots, 
i Pj ^ov an undamped vibrating system must be orthogonal, 

«•») 

The latent projectors for the simple matrix pencil Mp+K ■ 0 are ortho~ 

gonal since this latent‘>problem can be related to the eigen-problem. If 

2 

M ^ then 

K-My - )M^^^ (2.12) 

where it is assumed that M is positive definite. It then follows that the 
determination of the eigenvalues and eigenvectors of K-lp is a standard 
eigen-problem. The eigenprojectors of this algebraic system must satisfy 
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th« uiu«l orthogonal, piroparty with 
^Ki,0 ^KJ,0 “ ° 

Tha aigonprojnctora of K-Ih and tha latant projactora of K-fX^X ara ralatad 
with 


^Ki,0 




^ 21 - 1 ^ 21-1 ^ 21 - 1 ^ 21-1 


( 2 , 13 ) 


and 

only a linaarly indepandant latant vactoci for K-ly • tha latant vactors in 

A A 


Tha latent vactors of X^^. 


^21 equal and since thare ara 


Q and Q must ba orthogonal in ordar for (2,12) to hold with j^i. 
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3. The Cho leaky Decompoeitloh and the Canonic Fom 
Any aynooetclc positive matrix B m b can be decomposed into a product 
of two triangular matrices, Q and by the Cholesky algorithm. The canonic 
lambda matrix form can be obtained by factoring K out of (1*2), to the 
right or to the left, but this procees would destroy the syoBMtry of the 
resulting laid)da matrix. Bather than carry out this operation, the naea 
matrix can first be decomposed by the Cholesky algorithm [26^ with 


M- QQ^ 


(3.1) 


Ml HP 

where QeR is a lower triangular matrix and Q will be an upper triangular 
matrix. If Q is now factored to the left and to the right, (1.3) takes 
the form 


A(X) - Q[IX^+CX+K]q'^ - QA(X)q'^ 


(3.2) 


f l) -1~ -T -1 — -T ^ — 

where ' C ■ Q CQ and K Q 1C Q . Since C and K have been assumed 

to be symmetric, C and K will also be synmsitrlc. 

The latent roots of ACX) will he equal to those of A(X). which can be 

easily shown. Defining A(X) as. 


[A(X)r^- I[XI-Jr^ 


(3.3) 


where la the left latent vehtor matrix and is the right 


latent vector matrix, then 


I A(X)3 ^ 


(3.4) 


The superscript (-T) will be used to denote the transpose of the inverse 
of the matrix. 
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Th« noCTuilisation procatt gi\*n in C3*2) Srlll thnMforn mp th« InUnt vnc- 
fcosi of A(X) into with o 

- C(\ (3*5) 

In idditiont A(\) tuii ayMMitry thue 

- y^ ■ (3.6) 

Tho Choloiky dacomposltion and factorization will laava tha latant roota 
invariant but tha latent vectors will ba nodifiad as in (3.5) and (3% 6). 

The block companion matrix associated with A(X) will ba aa given in 
(1.5) and its aiganvtUuaa will ba equal to the latant roota of A(X). Tha 
right aiganvactor of A^ are given by 



and tha left aiganvector of A^ is 


(3.7) 


z 


cl " 





yi(Xii-W) 

^1 


(3*8) 


for with i - 1,3,5, ... ^m“l. Tha eigenvectors for the eigenvalue X^ with 
index i •• 2,4,6,. ..,m are equal to the complex conjugate of (3.7) and (3.8). 
It will be aasumad throughout the development that the latant roots of 
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A(X) afe diitlnct axccjipt: for thoaa at bhd origin^ X<^0» All latant 
roota will appaar In conjugata paira including thoaa at aaro which mat ba 
praaant in an avan nuariiac and all latant vactora will occur in conjugate 
paira. 

The i onpanion form aatrix haa the apectral reaolutlon 


(E.4) 


2u 

I ^ 

1-1 


1,0^1 


m 

I 

i-1 


^^i»0^l‘^^l,0^i^ 


and the partial fraction expanalon of (Xl-A) 


-1 


ia given by 


rA a)]“^ 

- c ' ■ “ 


[xi-A^r^ 



(3.9) 


(3.1Q) 


which can alao ba defined aa 


[A(X)r^-2 I ~ M i O_ R , k0 I fro i , 


1-1 X^+2a^X+(cjj4oJ) 


(3.11) 


where “ Re(Pj^ q) and q " ^ 0^‘ ^ (3*9) la the apectral re- 


aolutlon of A , then 
c 


A^O) - XI-A^ - I - 


(3.12) 


which can also be given aa 


m 


A, (A) - I 

° 1-1 


^ ^^i ,0'^’^i , 0^ ^^i.o^i^^i ,o^P ^ 


(3.13) 
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It follows froB (3.13) that 


1 E 


(E.l) <-h,oKo^ • 


and 


0 l”l 
-K -C J 


Racalling that an algenproj actor la daflnad aa 


T 

n H w w m 

^1,0 'ci -'cl 


hh 


t5ja.l+C) yl) 

" X i * 


for normallzad alganvactora , than y^ will ba proparly noraallzad and 


(3.14) 


(3.15) 


(3.16) 


1 

P * *■ 


^1^1 

f ■; 
U' 

l' : ■ 

It J 




[! 

It 



^1,0 

i 

■1 

A 

A . 

: 



^1,0^1 


whara P^ q la 

a latant projector of the 

f ’ 

(3.17) can alao ba darlvad from [A (X)] 

c 

i • ' ■ 


— 1 

XI+C I 

I - 

[A^(X)1 ^ 

^ - [A(X)1 ^ 


1 : 

c 

t 


-K Xl 


(3.17) 


(3.18) 


II 
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thui 




p\_o(Xii«> p\^o 


X^I+C 

-K 








(3.19) 


Equation (3.14) can be uaed with (3.17) or (3.19) to durtvo the reault 


(h 




(3.20) 


and from (2.15) 


1-1 " 


(3.21) 




(3.22) 


Equatlona (3.20) and (3.22) expresaaa the Invariances for the aaslgnisant of 
deeping In the ayatem and are key equations. The asalgnTRent of daaplng will 
be dl8ci,i9aed later. 

TWO of the properties of the elgenprojectors can be used to derive 
properties of the latent projectors. The elgenprojector are Idcnpotent 
matrices thus 
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Subitltutlng for ^ into thn nbovo aquation givnsi 


^.0 ■ - 'l.O ^.0 


(3.23) 


and 


(L.5) 


(3.24) 


Tha companion matrix will hava 2m linearly indapandant al^anvactora 
in which caaa tha aiganprojactora ara orthogonal with 

Vo'j.o-‘> 

Subatituting for ^ «nd P^ ^ and uaing C3.1'7) givaa 

(L.6) F^^q(X^1+1jI+C)Pj Q - 0 (3.25) 

and 

a,7) ° 

Properties of the eigenprojectprs for the companion matrix and the 
canonic lambda matrix are aumnarlsed in fable 3.1. These properties form 
the basis for the development of the assignment of damping. Each of thaaa 
properties have been verified with an example and a cotsputer rim. 
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" i 
t I 


5tl Siunary of fch* Propirf^^is of th« Elgan/lutont Projoctori 
(r dlitlnct rooti of multiplicity 


*• Block t!oirp«nion Form 


p^,i - mlgtnprojmctorij P^,A - latent projector 


=•1 l ' l-o ■ 1 
1-1 ‘ “ 


' l.o h,o ■ h.o 


ii " X|t2^«eiiir 


II 


^ i ,0 ^ j ,0 " ° 


1 ^ j 


E.4 P, . - (A - X.IiP- , ^ 
1 »J Cl 


j ■* Ifly • • tO^—l 


E.5 A_ 


-cj ■ Ji “’ l.O ' l,i 

1 n 

“I 


E.6 [A^( - [I X^ + CX + K] 


XI + C I 
XI 


I! 


E.7 P^^Q - {(X-X^) [A(X-X^) [A(X) ]”h| 


X"X ^^1 root) 

i 


E.8 P 


J -1 


i j " Jf” ~ P ' i t <■ ”l"^[A(X) ]“^>| ( X^ repeated root) 

d X i 


E,9 P 


1,0 


\,0 


“^ i . O *' 


E.IO tAjX)]*"^- I I 


m.-l 
r 1 


"i.O 

Vt.O 


i.£ 


i-1 A-0 ( X-X^ ) 


A +1 


li 


II 


r ^ 


I.J 


1 1 
5 ^ 
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* 

El XX 

Q ^ diwtlnct: tooc 

+ root 


• 

E.XZ 

*cAj ^ ‘ ^ f 



E.X3 

r 

A , ■ I A 

A i-1 =*• 



E44 

^cl “ ^1,0 




Aa) - IX^+CAfK Canonic Fora 



L.X 

2ni A 2 '** 

I ®*i 0 " ^ ® 

1-1 1-1 



Ii.2 

2m A m 

Jl ^1,0^1 ^^2i-'j.*0- ^21-i’^^2i-X»0^2i-X^ " ^ 



L*3 

2m A 9 *9 

j, <'i,o^i«i.o^i > ■ -<= 



L.4 

"l.O ■ 


‘ ^ 
i ^ 

X • 5 

■ 'l,0<^-«^1.0 



L.6 

- » 

1 ^ J 


L.7 

0 - 0 

1 f j 


L.8 

^ Ty. T . - -T 

^.0 ■ Vi/'!-! i\- J'li’ ^1^1 



i 


// 
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4* Spectral factorization of k(X) 

Saction 3 of thia chapter epnaidars tha raatrlctiona on tha aigan-- 
projactora and the latent projectora of tha coapanlon fora of A and the 

Cl 

canonic fora of A(l). In addition to thaaa raatrlctiona » atudlaa of A<X) 
will ravaal other raatrlctiona • 

Conaldar tha undaapad laabda aatrlx for ':ha atructura with 


. I 


A(A) 


IX^+K 


C4.1) 


i 


The aatrlx coafficienta of A(X) are aXiij and A(A) will have 2a latent roots 
Let K denote a mxa aatrlx with eigenvaluea where and 

aiallarly let K" have elganvaluaB X^ - with The fraa-fraa 

latent roota at X^^ ■ 0 will be equally apllt into the apactrua of and 
that of K • Tha laabda aatrlx adalta the apactral factorization 


I i 


A(X) - (IX+kT) (IX+K“) 


C4.2) 


ii 


where for the undaiqpad caae of (4.1) 


K - k'^k" 

0 - k*+k" 


(^.3) 

(4.4) 


The atlffnaaa aatrlx haa alganprOjactora thua K can be conatructed by 

the aua 


g t 

H 


a 


‘‘ki.oV 


».5) 


I J 
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n X 

wh*rt ■* Nov K can b« factortd^ «• glvtn in (Ai3) with 
nnd 


^ ^Ki ^K1 " ”^*^1 

4* « 2 

vlth » X|y^ •■ o)^t Notn thnt (4.4) !• intlailad ■lnc« 

K+«- - tfK1.0J“r"KI,0J“i) - 0 


(4.7) 


(4.S) 


aa raqulrad. An Bxaaq>la can ba glvan to illuatira^:;a tha fac/iorizatlon 



K(I»J) MATRIX 


9 

^'5 

0 

-5 

11 

-6 

0 

—6 

13 

EIGENVALUES 

3.11262 

10.6609 

19.2264 

EIGENVECTOR MATRIX 

.587523 

.737925 

.332087 

.691794 

-.24513 

-.679214 

,419805 

-.62879 

.654513 

EIGENPROJECTOR FOR LAMBDA - 3.11262 
.345184 .406446 

. 246645 

,406446 

.47858 

.290418 

.246645 

.290418 

.176236 

EIGENPROJECTOR FOR LAMBDA - 10.6609 
.544534 -.180888 

-.464 

-.180888 

.0600886 

. 395377 

-.464 

.154135 

.395377 


t Note that -jlC*' and jx“ rta aquara roota of K, l.e. K^'^. 


i 
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EIUENPROJECTOR iFOR LAMBDA - 19.22( 4 


.110282 

-.225558 

.217355 


-.225558 

.461331 

-.444554 


.217355 

-.444554 

.428387 



Using (4.6), Is given by 



2.87052 

-0.862567 

-0.12.i807 


-0.862567 

3.06338 

-0.933635 


-0.126807 

-0.933635 

3.48026 


end k“ - -K*". 


When dsiaplng Is to be added, K Is Invariant under damping assignment 


thus the lambda matrix becomes 


Ajj(X) - IX^+CX+K - (IX+kJ) (IX+kJ) 

where 


C4.9) 




C4.10) 



C4.ll) 


The Invariance of K can be satisfied by postmult Ip lying by a unitary 
matrix and premultiplying K~ by the transpose of the conjugate of the uni- 
tary matrix, l.e. 


K 


+ * - 
K UU K 



C4.12) 


from ^>hlch It follows that 
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c » K^U+uV » k‘*‘u+w*(k‘^)* 


(4.13) 


Th« laost general form of U la complex with U ■ *lnce la 

Imaginary with K* « ~K^ 


C - (K^+j 4)(Ug^-j Uj)+(uJ+J Ui)<kJ;-J 4) 


(4.14) 


where " 0 for the undamped caae therefore 


Re(C) - 0 - K^Uj+UjK^ 


Im(C) - 0 - KjUj^-uJKj. 


(4.15) 

(4.16) 




Recalling that K » 


K - kJk^ - CUj^-j Uj) (Uj^+j Uj) V 


and that K - (jK^^^) (-jK^^^) - kV, (4.17) takes the form 


(4.17) 


K - K^^^(U^+j Uj^)V''^ " 4' S 


(4.18) 


The upper-half spectral matrix must have eigenvalues equal ti) the latent 
roots of A(X) for the upper-half plane. Similarly, must contain eigen- 
values only in the lower-half plane. Figure 4.1 below shows the eigenvaliie 
locations for the undamped and damped matrices, and K^, 


ri 
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eigenvalues of damped 
system 








X* 


X.. 




■*5k 


eigenvalues of undamped 
system 


Figure A.l 

The damping problem requires that a unitary matrix be found such that 


U - Uj 

U(U*) - I 
T T 

Vr '’i“i - 1 

T T 

Vr “ Vi - ° 


(4.19) 

(4.20) 


(4.21) 


(4.22) 


where U* is the complex conjugate transpose of U. The das^tlng matrix is 
given by 


C - Uj, + Uj 


(4.23) 


m 


with th« additional conatralnt that C ba triadiagonal and that 


0 




(4.24) 


It follows from the above analyais that the lambda matrix A(l) admits 
a spectral factorization. 


A(X) - (IX+k'^)(IX-W(“) (4.25) 

for the undamped case with K ■ and 

Ajj(X) - [IX+kVj^-J U3.)][IX-KuJ+J U^)k’] (4.26) 

for the damped case. The complex unitary matrix must have the property that 
C as defined in (4.23) properly accounts for the damping liqilemented in the 

nXffi 

structure. If C is trldlagonal then UeC must satisfy the tridiagonal 
proparty as well as satisfy (4.24). 

The matrix function 

C - f (U) - (k’^)U + U*(JC*‘)* 


given in C4.13) has the following properties; Taussky and Wielandt [28] 

a) The nxn complex matrices and U are nonsingular if f(U) Is a 
positive definite matrix 

b) f(U) is linear over the subfield R of real numbers. 
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Property a) con be shown «• follows; 1st z dsnots s vsctor such that Uz ■ 0 
which Implies that U is singular. It follows that 

K^Uz-fz*U»>CK.‘*’)''' - **fCU)z - 0 (4.27) 

4. 

Thus f(U) is not positive definite. Similarly 1 if K is singular then 
zK^ ■ 0 and 

ZK.‘^U+U*CK^)*Z* « z fCU) z* - 0 
and this Is a contr^)idlctlon of a) . 

Property a) then Implies that if K is singular, e.g. eigenvalues at 

the origin then C will be singular. Since K con be positive seml-definite 

4. 

the C will be positive semi-definite whenever K is semi-definite . The 

4. 

Buitrix K is the spectral factor of K thus K will also he positive seml- 
deflnite when is positive semi-definite. 


chapteI 6 

OPTIMAL CONTROL OF SELECTED MODES ,, 

;Y 

The usual approach taken In applying optimal control theory Is to 
determine state or output feedback for all of the modes of the system. 

This procedure works quite well when the number of modes in the system is 
not large but the computational load for several hundred modes makes this 
type of control impractical if time varying gain is used. Even when a 
constant gain liV Used, ch.e computation of the gain is not a trivial task. 

There have been numerous papers, see [29]- [30], published in applying 
optimal control theory to structure with the development based on reduced- 
order models. The computational load can be reduced significantly by this 
approach but the reduced-order model must be carefully chosen if mode spill 
over is to be avoided. 

The work in this section will take an entirely different direction. 

The computational load for the procedure is reasonable and the mode spill 
over problem can be eliminated. The spectral factorization algorithm will 
be used to decouple the selected modes from other modes of the structure. 
The optimal control theory will then be used to construct the feedback for 
the selected modes. The uncontrolled modes are uncoupled from the control 
modes and the possibility of mode spill over is eliminated. 

The method presented in this chapter is similar to that of Oz and 
Meirovitch [31] where the modes are decoupled and the bptimal control is 
determined for each mode. The major difference is in the method of de- 
coupling and the computational procedure. 


113 




-^-4^ -^=^'-'4 


1, Introduction 

Consider the state equation for a plant with 


XU 




X(t) - AX(t) + BU(t) 


( 1 . 1 ) 


nVn v%^ffti 

where AER” ” is the system matrix, BeR is the input force matrix, 
X(t)eR”^^ is the state vector and U(t)eR™^^ is the force vector. Let J be 
the associated scalar cost function with 

J(X,U,t) - ~ x'^(0)HX(0) -f “ [x'‘^(t)q^X(t)4‘U^(t)Q2U(t))dt . (1.2) 

The Hamiltonian for the system is 

M(X,U,P,t) - i x'^(t)Qj^X(t) -f I u’^(t)Q2u(t)+P^(t)[AX(t)+BU(t)) (1.3) 

from which it follows that X(t), P(t) and il(t) must satisfy the equations 


X(t) - AX(t) + BU(t) 
P(t) - -Qj^X(t)-A'^P(t) 

0 - Q2U(t)+B'’^P(t) . 


(1.4> 

(1.5) 

( 1 . 6 ) 


The desired control for minimizing the cost function is 


U(t) - -Q2Vp(t) 


(1.7) 


where it will be assumed that 


P(t) - RX(t) . (1.8) 

Differentiating (1.8) with respect to t and using (1.4) anii (1.5) gives 
the algebraic Rlccati equation 

+ a"^R + RA - RBQ2^B^R - 0 


(1.9) 


// 
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for R. Subotltutliig R into (1.8) and tho resulting equation into (1.7) 
gives the control 

U(t) - -Q2VRX(t) (1.10) 

This control vrlll give the closed loop matrix 

I - A-BQ^^b’^R Cl. 11) 

It is usually assumed that H, and are symmetric positive semi-definite 
matrices and Qg symmetric hut positive definite. The Riccati matrix 
obtained from (1,9) will also be symmetric and positive definite. The 
matrices H, and Q 2 are weighting matrix chosen to fix the cost penalty 
for tbe initial conditions , the displacements and the control effort. 

An example will now be given to illustrate the computational procedure. 
Let A be 


0 0 10 



H - 0 

■Jl - ^4X4 
^2 " ^2x2 


Substituting into (1.9) gives 


9 
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10 0 0 


0 0-9 r 


0 ( 0\ 1 0 

0 10 0 
0 0 10 
0 0 0 1 

+ 

00 5-11 

10-10 
0 10-:^ 

R+R 

Q 0 0 1 

-9 5-10 

_5 -11 0 -1_, 


0 0 0 0 
0 0 0 0 
0 0 10 
0 0 0 1 


The algebraic solution to this equation Is 


R - 


4*1817 

-2.02521 

0.0736876 

0.033291 


-2.02521 

4.99179 

0.033291 

0.060371 


0.0736876 

0.033291 

0,465215 

0.022792 


0.033291 

0.060371 

0.022792 

0.456098 


and the closed loop matrix Is 


A ■ 


0 

0 

0 

0 

0 

0 

0 

1 

-9.07369 

4.96671 

-1.46522 

-0.022792 

4.96671 

-11.0604 

-0.0227921 

-1.4561 
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The eigenvalues for the open loop system are X - -0. 51 j 2, 15661 and 
X -0.51J3.85344 whereas the closed loop eigenvalues are X ■ -0.7419521 
j 2. 10985 and X - -0. 7187081 j 3. 82303. 

It should be noted that If the closed loop system matrix represents a 
model of the form 


A ■ 


0 

-K 
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It It ntcesttry to modify the ttlffnege nmtvlx;. It should tlto be pointed 
out thtt G does not represent e model with pestlve demping es C does not 
heve the proper structure. If the closed loop system hss the fsedbeck de- 
fined In (1.10), there are no restrictions since this control lev is not 
for a passive system. / 

The control vector for optimal control should I always be constructed 
from the velocity elements of the structure if the matrix form given above 
is to result from the control! Basically this means that the B matrix 
should have zero elements in the upper half, i.e. in the first nXm block 
of B. 

The algebraic Rlccati solution for the exanqp>le given above required 
that a 2nx2^; matrix be used where A is nXn. It is obvious that a system 
with a large number of modes is not suitable for the above type of analysis. 
The optimal control procedure above will also add damping to all modes 
where it may only be necessary to damp e few modes. 

The analysis presented above is more difficult when an undwnped system 
is considered. The reason for this is due to the presence of multiple 
eigenvalues of the matrix considered in solving the algebraic Rlccati 
equation. The mode decoupling method presented in the next section will 
consider an undamped system and the cooqmtatlonal problems with the alge- 
braic Rlccati equation will be eliminated. 
i 


lU 


2, Ho<to Decoupling 

A method will be given in this section that decouples some of the 
nodes from the remaining ones so that the optimal control atrateg/ can be 
carried Out on a lower-order system. 

f 

Consider the undamped system matrix with C’^O such that A is 



( 2 . 1 ) 


where K is positive definite. The eigenval\^es of A are along the Ju) axis 

and occur in complex conjugate pairs. The eigenvalues of K are given by 
2 

where is an eigenvalue of A. This suggests that the spectral deconr* 
position of A can be obtained from considering K rather than A. 

Suppose th.e there exiet. . elmiletlty trenstorBetioe metrlx auch 

that 


Vk" ■ 


Si « 

« Si 


( 2 . 2 ) 


Where has eigenvalues X^|p|and Kg 2 has eigenvalues Xj^>|p|. If such a 
matrix exists, then gives the spectral decomposition of K. To find T^^, 
the eigenvectors of K must be found or the sign algorithm can be used to 
generate T. The eigenvector procedure will probably be the most efficient 
for large systems so the procedure will be based on that method. 

Let denote the eigenvector of K from which it follows that the in- 
verse of T is Constructed from the columns of Consider the stiffness 
matrix 


K - 


9 -5 

-5 11 

0 -6 


0 

-6 

13 


1X0 
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for which the elgenveluei ere 


Xj_ « 3.11262 
X 2 - 10.6609 
X 3 - 19.2264 


and the eigenvector matrix ie 




0,587523 

0,691794 

0.419805 


0.737925 

-0.24513 

-0.62879 


0.332087 

-0.679214 

0.654513 


The eigenvectors are ordered in the same order as th® i'tgenvelues . 

Assume now that is to be constructed such that contuins Che 
eigenvalue X^^ and K32 eigenvalues X2 and Xg. The eigenvector matrix 
is partitioned with 


'*'11 

mmi0i 

*12 


"0.587523 

0.737925 

0.332087 

- 

0.691794 

-Ov 24513 

-0.679214 

't'21 

'*’22 


0.419805 

-0.62879 

0.654513 


-1 

and T is then given by 


2 


-I 


"'^21‘^li 




or numerically 


-1 

12'^22 
I 


.-1 


-0,5 -0.58874 -0.357266 

-0.59874 0.5 0 

-0.357266 0 0,5 


The Inverse of Tj^ is 


<2.3) 
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i: 



- 0.690368 

-0.812892 

-0.49329 


0.812892 

1.04284 

0.580837 


-0.49329 

-0,580837 

1.64753 
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The almllarlty tranefom as expresaed in (2.2) gives 



3.11262 

0 

0 


■S.1 

0 


0 

16.8874 

-2.42734 

M 

V 

0 

-6 

13 


0 



which has the correct eigenvalues. 

If is to be 2x2 yith eigenvalues and X^. the partitioning of 
would be changed so that is 2X2, ^22 transforma- 

tion matrix is then constructed and the similarity transformation applied. 

Thus far, the spectral decomposition of K has been carried out but 
the system matrix given in (2.1) must be considered as this is the matrix 
that is of concern. Let T be a new transformation matrix with 

(2,4) 

Using T as given in (2.4) then 
TAT"^ » 

gives a new system matrix with in the lower loft corner of the matrix. 

Substituting the value of from the example, the new matrix Is 




0 T, 


K 


0 I 
-K 0 




0 T, 


K 


-1 


-Tj^KT"^ 0 


(2.5) 


T 


\ 0 




0 
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which l0 not block dlagonaXlased. To block dlagonaliza TAT**^, conaCruct a 
£^ow- column Interchange matrix E where 



The block diagonal form can then be found with 
■ ETAT“^E - T^AT^^ 



The spectral decomposition process will modify the state vector X(t) 
as will now be shown. Let V(t) be defined as the transformed vector 


V(t) - T^X(t) (2.7) 
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therefore 

V(t) - T^XCt) - T^AX(t) . (2.8) 

but X(t) - T"H(t) thu3 

V(t) - T^AT^^VCt) (2.9) 

The similarity transformation on A to block diagonalize A will map X(t) into 
a new vector V(t) where 

C2.10) 


( 2 . 11 ) 

ET* is 
A 

0 0 

-0.690368 -0.812892 

0 0 

0 0 

-0.812892 1.04284 

-0.49329 -0.580837 


v^(t) - -0.690368x^(t)-0.812892x2(t)“0.49329x3(t) 
V2(t) - -0.690368xj^(t)-0.812892i2(t)-0.49329x3(t) 

Vg(t) - -0.49329 x^(t)-0.580837x2(t)+1.64753x3(t) 


P 

-0.49329 

0 

0 

-0.580837 

1.64753 


V(t) - E 


‘K 


X(t) 


K 


or 


XCt) - 


[.”1 

0 


E V(.t) . 


T, 


K 


For the exanqile given in this section T 


T - ET, 


-0.690368 

0 

-0.812892 
«0. 49329 
0 
0 


-0.812892 

0 

1.04284 

-0.580837 

0 

0 


-0.49329 

0 

-0.580837 

1.64753 

0 

0 


thus the vector V(t) has components "'^.(t) with 
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T • • • 

where X (t) »• [Xj^(n) X2<t) X3<t) Xj^(t) X2(t) x,^(t)] for the undamped model. 

All of the computations for the decomposition given In this section 
are carried out by considering the K matrix which is mXm, It is not neces- 
sary to find the eigenvectors of the A matrix since the necessary Information 
is contained in K. ‘ 



3. Optimal Control o'i tK:^ Undamped Decoupled State Matrix 
It vas shown In the previous section that the state matrix could be 
block dlagonalliijed with selected eigenvalues o£ A placed In one o£ the 
selected block matrices, let the block matrix £or the undamped syatem 
have the general fotm 


( j' 

(3.1) 

where has eigenvalues |X^|<P and A^2 has eigenvalues with p a 

scalar variable and A Is the undamped matrix. The value o£ p will be 
chosen to Include the desired modes In A^^> 

Consider now the algebraic Rlccatl equation lor R and let A^ be the 
decoupled matrix, thus R must satls£y 

Q^+aJr+RA^-R B q"^ R - 0 (3.2) 



where and are weighting matrices £or V(t) and U(t). The matrix B 
represents the control Input matrix where 


V(t) - A^V(t) + BU (3.3) 

with Ag de£lned In (3.1) and B TB. It will be assumed that the algebraic 
Rlccatl equation Is completely decoupled such th^t 

*^12 ^2^2 *^^2 ” *^2^® *^2 ® ^2 \ ® (3.5) 


The Rlccatl matrices R^ and R 2 can be £ound Independently since the equations 
are decoupled. 
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Substituting for In (3,2), 

\\ 

+ 'fWj + rtat“^ - r b r - 0 ( 3 * 6 ) 

! 

and rearranging gives 

+ a^t'^t + t’^ta - t^rtbq“^b'^t^rt - 0 (3.7) 

Defining R ■ T^T and gives 

+ A^R + RA - RBQ“^B^R - 0 (3.8) 

which is the usual algebraic Rlccatl equation for the general optimal con- 
trol problem. 

Denoting equation (3.4) as system 1 and (3.5) as system 2, It follows 
that system 1 has the system equation 

Vi(t) - A3^V^(t) + B^Uj^(t) (3.9) 


with cost 




r°° 

[V^(t)Q^^V^(t) + u^(t)Q^,u,(t)]dt 


*2rv 


(3.10) 


The other System has the state equation 


V£(t) - Ag2V2(t) + B2U2(t) 


(3.11) 


and cost 


^2 ^^2 ’“ 2 *^^ 


[V2(t)Qj^2^2^^) + U 2 (t) 022^2 


(3.12) 


where the Initial state cost has been neglected. Assuming that the first 
system Is the desired system for damping, then B 2 " 


0 Will leave system 2 
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undamped R 2 ■■ 0. It then follows that the uncoupled system Rlccatl equation 
is 


R - T* R T 


^lA'^ll ’^llVl2 


^12Vll ^^12^12 


(3.13) 


Since Rj^ will be symmetric then R will bo symmetric as desired. The control 
input vector ¥ will have the form 



1 

I 


h~ 



m 





0 


(3.14) 


therefore 


B - 


-I— 

T B - 


(Tii-T^^r^Jr^^) ^^21~^22‘^12'^11^ 


11 *12"22 2r 

^12"^ll'^21^22^ 


(T, ,-T, ,t;Jt„)“^ <^22“^2 A^12^"^ 


(3.15) 


The numerical value of B^^ can be chosen such that B^ is zero and B 2 y* 0. 
The closed loop system matrix is then given by 


A - A-BQ2^b\ 


0 I 


0 

0 

-K 0 


^2®2'^12\^11 

®2»M2Vi2 


(3.16) 


,T„T 


where it may be possible to make B2B2'^12Vll " ° properly selecting B 2 
and the weighting matrices. In general, this matrix will not be zero and 
the stiffness of the structure will be changed. 

The example used in Section 2 will b(fe used to illustrate the computa- 
tional procedure. The stiffness matrix was 


K 


9 -5 0 

-5 11 -6 

0 -6 13 


I ' i 


f I 


: ! 
! 5 


i i 
i , 




I ; 


i i 




n 

t > . 
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wheve 



and 


0 1 
j-3. 11262 


0 

0 

1 

0 

0 

0 

0 

1 

-16.8874 

2.42734 

0 

0 

6 

-13 

0 

0 






The flrat step is to select B such that 0 and B^ “ 0 which can by ob- 

tained from fl - TB where 


-0.690368 

-0.812892 

-0.49329 

0 

0 

0 


^11 

0 

0 

0 

- .690368 

- .812892 

-0.49329 



-0.812892 

1.04284 

-0.580837 

0 

0 

0 


^13 

-0.49329 

-0.580837 

1.64753 

0 

0 

0 


^4 

0 

0 

0 

-0.812892 

1.04284 

-0.580837 


^4 

0 

0 

0 

-0.49329 

-0.580837 

1.64753 


h5 


0 

1 



L“J 

The solution to this equation is 

B^ - [0 0 0 -0.5 -0.588738 -0.357266] 

-1 

where the last m elements is the vector of belonging to the controlled 
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node* 


Selecting the 2x2 metrix for damping ^ then for ■ I, ■ I 
— T 

and •* [0 1], the Rlccatl equation Is 



which has a solution 



and Agj closed loop is 


■■ mrnt. MM a. 1 

Si " Sr®i‘^2 W " 


0 1 
-3.2136 -1.09363 

mmm 


recover the 

full Rlccatl matrix, 

use R 

T^fe which 

gives 

2.61365 

3.07751 

1,86753 

.0481279 

.0566694 

.0343889 

3.07751 

3.62369 

2.19898 

.0566694 

.0667269 

.0404921 

1.86753 

2,19898 

1.33441 

.0343889 

.0404921 

.024572 

.0481279 

.0566694 

.0343889 

.521233 

.613739 

.372438 

.0566694 

.0667269 

.0404921 

.613739 

.722664 

.438536 

.0343889 

.0404921 

.024572 

.372438 

.438536 

.266119 


— -IT 

The closed loop system matrix is given by A ■ A-BQ^ B R- where the numerical 
values are 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

-9.034S6 

4.95896 

-.0249062 

-.377504 

-.444502 

-.269739 

4.95896 

-11.0483 

5.97067 

-.444502 

-.52339 

-.317611 

-.0249063 

5.97067 

-13.0178 

-.269739 

-.317611 

-.192737 
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The eigenvalues of A are 

h 2 '“ -0 • 5468161 j 1.70722 

- ±j3, 26511 

- ±j4.3848 

which agrees with those of and Aq 2 » 

The required feedback control vector u(t) Is given by 

UCt)--Q2^B\x(t) - -Q2 ^B^’r T X(t) (3.16) 

which can now be determined as Q 2 » ^ known. 
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4. Opfcloal Control of tho General Canonical State Matrix Fom 
The procedure described in the proceeding sections considered the un- 
damped state matrix in canonical form. The method can be extended to the 
general form by utilizing the algebraic theory in Chapters 2 and 5. Con- 
sider the damped matrix 



(4.1) 


where and CeR*“^® with no restrictions placed on the forma of these 

matrices. The optimal control theory can be applied to this form with some 
modifications in the computational procedure. 

The right and left eigenvectors for the eigenvalue are giv€«i by 





(4.2) 


and 


'cj 




(4.3) 


T 

respectively where 4(Xj)y^j ■ 0 and A (X^)z^j ■ 0, The eigenprojector 


for the eigenvector X^ is defined as 


P . , 7 zT 

j T ^cj cj 

j z . y . 

cj -^cj 


(4.4) 


where y^j , and z^^ denote the unnormalized eigenvectors and the eigenvectors 
with the bar indicates normalized eigenvectors. The eigenvalues will occur 



X3l 
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in complex conjugate pairs thus ^2i-l taken as the eigenvalue in 

the upper half plane with and lower half plane « 

The eigenp rejector for a pair of conjugate eigenvalues of the jth mode 

^j “ ^2i‘-l‘*’^2i " yc21-1^2i-l'*^c2i^ 21 

« 

or from (4.2) and (4.3) 


^21-1^2 i'-1^^2i-l^'^‘^ ^’*■^21^21 ^^21^’^^ ^ 

y2i-iyL-1^^2i-l^'^^ ^2i-l^'^2iy21^^21^’^ ^2i^ 


^2i~l ^21-1 ^21 ^2i 

y 2 i-l^L- 1 ^ 2 i-l’^ 2 i^ 2 i^ 2 i 


" ^2l-l'^^2i-l 


(4.6) 


but X 2 i - Xji^i and Thus the terms in (4.6) can be combined 

if desired. The elgenprojector will be a real matrix since. Aj ■ AF^ 
will be a real matrix. 

Su;:’P 08 e now that the first of modes, j " l,2,...,q, are to be the 
ones that are to be damped. It then follows that 


q 

^d’ j: 

^ 1-1 J 


and if is the undamped projector, it follows that 


(4.7) 



(4.8) 


The projector P^ can be conq>uted by determining the latent vectors of 
2 

A(X) - IX +CX+K or from the eigenvectors of A. 

The matrix T - S-fJ is needed but this can be computed from P^^. The 
elgenprojector for the first q modes can also be defined in terms of the 
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in 


algti futrlx S with 


Pjj - I (S4-I) - ♦ 


I 0 

q 


with 


S » $ 


0 -I 


m-q 






(4,9) 


(4,10) 


where $ is the eigenvector nuitrix of A, To determine S+J, note that 


S - 2Pp-I 


C4.U) 


and thus 


T - S+J - 2Pjj-I+J - 2Pjj+2 


0 0 

0 -1 

L. n-y 


(4.12) 


The next step in the algorithm is to find which is the block dia- 
gonal matrix that is to be damped by the optimal control strategy. It can 
be shown that 


" ^11 ^ 12^21 


(4.13) 


where and Aj ^2 partitioned blocks of A with and A^ 2 ^^^ 

These submatrices are determined directly from A, The Rlccati function R 21 
must be found but this matrix can be obtained from Pjj or T, It can be 


shown that R 2 ^^ is given by 


.-1 


^21 “ ‘^ 21*^11 

where and Partitioning (4.9) and carrying out the 

algebraic steps gives 


(4.14) 
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sign matrix S with 


Pjj - f (S+I) - 4> 


I 0 
0 0 


(4.9) 


with 


S - 4 


I 0 

q 

0 -I 




(4.9) 


m-q 


where 4> Is the eigenvector matrix of A. To determine S-fJ, note that 


S - 2Pjj-I 


(4,10) 


and thus 


T - S+J - 2Pjj~I+J - 2Pjj+2 


0 0 
0 -I 


(4.11) 


m-q 


The next step In the algorithm Is to find which Is the block dia- 
gonal matrix that Is to be damped by the optimal control strategy. It can 
be shown that 

^1 " ^11 ^12^12 • (4.12) 

where A^^^ and partitioned blocks of A with and 

These submatrices are determined directly from A. The Rlccatl function R 2 j^ 

must be found but 
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'd ■ 


♦aiWxr*12*M*21>'' *21»21-*22+«*U>'' 


(4*35) 


*21 • *B21 *DU 


C4.16) 


and 


^12 * “^Dll ^D12 * 

The latter siatirlx, Rj_ 2» is not needed unless T 
-1 

T can be written down directly with 


C4.17) 

Is sought In which case 



(4*18) 


The computational procedure for the damped system matrix is as follow: 


a) Select the modes to be damped* 

b) Compute the right eigenvectors for the selected modes. 

c) Using the latent vectors obtained from the eigenvectors^ conqpute 
the elgenprojector for the selected modes* 

d) Compute T and R 2 j^* 

e) Form the block matrix 

f) Using the laub algorithm [32] compute the algebraic Rlccatl solu- 
tion for the control law of 

g) Compute u(t) « -Q^^^RX(t) and Implement the control. 


It Is assumed that la not singular. 
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5» Gonclufion ^ 

The developnent o£ « procedure for optinel control of eelected nodei 
of e second-ord'^r lembde netrlx has been given in thia section. Although 
the development has been abbreviated and did not cover the most general 

4. 

model there appears to be no significant problems for the general case. It 
has been shovm that the control vector for a rather amall aacond-order 
lambda matrix can be determined so that damping is added to the lowest 
modes. The other modes remain undamped. 

The mode decoupling procedure used with the optimal control strategy 
avoids the problem of mode "spill-over". The algorithm as presented does 
not give the design the option of placing system eigenvalues at desired 
location although this can be achieved with some modifications to the theory 
provided that the designer give up some restrictions on the final form of 
the closed loop system matrix. 

The calculations carried out In this chapter resulted from several 
specific computer programs which could be combined for a general package. 
Since software development is a time-consuming task» no effort was made 
the develop such a package at this time. 


ORIGINAL PAGE B 
OF POOR QUAUTY 


136 


Appendix A 


The elgenprojectors ^ (As defined in (3«2)(Can be obtain- 


ed from 


^in 


l<k) 

dl* 


^i,m^-0.^-k " ^ (AI-J)”-"1}W^ (A.l) 


®i“*li'^^ -1 1 T 

k * a. /I T T\ +n\»jA 

L 


k " Of 1 f t • • (iB^— 


where W» and W, are the right and left eigenvector matrices. The inverse 
K L 


of (Xl-J) is given by 


(M-J) 


-1 




-1 


(XI-J2) 


-1 


(XI-J,.) 


-1 


(A. 2) 


m±Xm. 

where J^cC ^ is the Jordan block associated with the eigenvalue X^ of 
multiplicity m^. For a defective matrix, the Jordan block J^ can have the 
following structure 


" 1 - 


•\ 

I 

' I 

,> 

I 
I 




qi-1 


X., 1 / 

■i 
I 
I 




1 


(A.3) 


mj-qj+l 
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whsre is the nuDber of prJiunry elgenivoctori and the number of gen- 
eralized eigenyactora aaBOclated with X^. It followa that la 

given by 


(XI-J^) 


-1 




-1 




-1 


(X-\> 


-1 


<X-X^)’^ (X-X^)'^ ... (X-Xj)’^'"^'* 

,-l 


(X-X,) 


L 


The bracketed term of (A.l) for (Xl-Jj^^) gives 


(x-xp 


-1 


(A. A) 


®4“q4'^x 

(X-X^) ^ ^ (XI-J^) ^ 


(X-X^) 


"r^i 


(X-X^) 


«l-qi 


(X-X^) 


V**! 


• (X-X,) 




(A.5) 


It la not difficult to evaluate (A.l) for k ■ 0,1» . . . »m^-q^. The elgenpro- 
jector ^ Is then given by 




^»V»1 " \ 


u 

0 


V “1] 

o! 


1 * - 
I 0 


■ “r •" '>l 


(A. 6) 


.1 


w 

X38 
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The shifting matrices are important in computing eigenprojectora of a 
matrix provided that the eigenvalues and eigenvectors are known along with 
the multiplicity of each eigenvalue and the degeneracy of each Jordan Block. 
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APPENDIX B 

The Kronecker product is a useful mathematical tool for aystema theory 
and can be defined as follows. Let and then the Kronecker product 
A©B Is 


A©B - 


An excellent treatment of Kronecker algebra has been given by Brewer [14]. 
If the matrices 0^^^ and are considered the following pro- 
perties are of Interest In this paper: 

1) (N©M)"^ - n”^©m"^ 


“ll® “l2® ' ’ “iq® 


* 21 ® 


a -B 

pi 


a B 

pq 


(B.1) 


11) (A©B)(C©D) - (AC)©(BD) . 


<B.2) 
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